arXiv:1502.06210v2 [math.DS] 24Jul2015 


REGULARIZATIONS OF TWO-FOLD BIFURCATIONS IN PLANAR PIECEWISE 

SMOOTH SYSTEMS USING BLOWUP 

K. ULDALL KRISTIANSEN AND S. J. HOGAN* 


Abstract. We use blowup to study the regularization of codimension one two-fold singularities in planar piecewise smooth 
(PWS) dynamical systems. We focus on singular canards, pseudo-equlibria and limit cycles that can occur in the PWS system. 
Using the regularization of Sotomayor and Teixeira m, we show rigorously how singular canards can persist and how the 
bifurcation of pseudo-equilibria is related to bifurcations of equilibria in the regularized system. We also show that PWS limit 
cycles are connected to Hopf bifurcations of the regularization. In addition, we show how regularization can create another 
type of limit cycle that does not appear to be present in the original PWS system. For both types of limit cycle, we show 
that the criticality of the Hopf bifurcation that gives rise to periodic orbits is strongly dependent on the precise form of the 
regularization. Finally, we analyse the limit cycles as locally unique families of periodic orbits of the regularization and connect 
them, when possible, to limit cycles of the PWS system. We illustrate our analysis with numerical simulations and show how 
the regularized system can undergo a canard explosion phenomenon. 

Key words. Piecewise smooth systems, blowup, geometric singular perturbation theory, sliding bifurcations, canards, 
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1. Introduction. Piecewise smooth (PWS) dynamical systems 15, [27j are of great significance in 
applications [8], ranging from problems in mechanics (friction, impact) and biology (genetic regulatory 
networks) to control engineering [32]. But, compared to smooth systems [16], the study of PWS systems is 
in its infancy. For example, notions of solution, trajectory, separatrix, topological equivalence and bifurcation, 
all need revision and extension m ■ Often PWS systems are used as caricatures of smooth systems Ei Eg, 
especially if significant amounts of computation are expected. So one of the major challenges of PWS system 
theory is to see just how close the behaviour of a PWS system is to a suitable smooth system. 

In this paper, we focus on PWS systems in the plane, of the form: 

x = A ± (x), x e E± C R 2 , (1.1) 

where the smooth vector fields X ± , defined on disjoint open regions E±, are smoothly extendable to their 
common boundary E. The line E is called the switching manifold or switching boundary. The union E U 
E_ U E+ covers the whole state space. When the normal components of the vector fields on either side of 
E are in opposition, a vector field needs to be defined on E. The precise choice is not unique and crucially 
depends on the nature of the problem under consideration. We adopt the widely-used Filippov convention 
H9> where a sliding vector field is defined on E. In this case, the dynamics is described as sliding and the 
PWS system m, together with the sliding vector field, constitute a Filippov system. Such systems possess 
many phenomena that are not present in smooth systems; grazing and sliding bifurcations, period adding 
bifurcations and chattering are (almost) ubiquitous in and (virtually) unique to PWS systems. 

Sotomayor and Teixeira m proposed a regularization of a planar PWS dynamical system, in which 
the switching manifold E is replaced by a boundary layer of width 2e. Outside the boundary layer, the 
regularization agrees exactly with the PWS vector fields. Inside the boundary layer, a monotonic function 
is chosen such that the regularization is at least continuous everywhere. The regualization of PWS systems 
in R 3 was considered by [26] and in R" by [24] , 

It is natural to ask whether bifurcations in PWS systems are close to bifurcations in a suitable smooth 
system. But for any regularization, there is a fundamental difficulty when dealing with bifurcations. Fenichel 
theory mmmm, the main tool used to analyze regularization, requires hyperbolicity, which is lost at 
a PWS bifurcation. A widely used approach to deal with this loss of hyperbolicity is the blowup method, 
originally due to Dumortier and Roussarie HI Holm, and subsequently developed by Krupa and Szmolyan 
m in the context of slow-fast systems. 
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Buzzi et al. [2] considered how different PWS phenomena]]] in the plane were affected by the regularization 
method of Sotomayor and Teixeira |3U| . A similar study in R 3 was carried out by Llibre et al. fU'll . 
Regularization of PWS systems in R” was considered by Llibre et al. [21]. These three papers considered 
the case of one switching manifold separating two different smooth vector fields. Regularization in the case 
of two intersecting switching manifolds was considered by Llibre et al. [25], and in the case of surfaces of 
algebraic variety by Buzzi et al. [3]. Regularization of codimension one bifurcations in planar PWS systems 
was considered by De Carvalho and Tonon [6]. Common to all of these studies, however, is that they do not 
deal rigorously with the loss of hyperbolicity at a PWS bifurcation and hence they do not properly unfold 
the effect of the regularization. 

Recently, Kristiansen and Hogan jTBj successfully applied the blowup method of Krupa and Szmolyan 
m to study the regularization of both fold and two-fold singularities of PWS dynamical systems in R 3 . For 
two-fold singularities, they showed that the regularized system only fully retains the features of the PWS 
singular canards when the sliding region does not include a full sector of singular canards. In particular, 
they showed that every locally unique singular canard persists the regularizing perturbation. For the case of 
a sector of singular canards, they showed that the regularized system contains a primary canard, provided a 
certain non-resonance condition holds and they provided numerical evidence for the existence of secondary 
canards near resonance. Other authors [29] have used asymptotic methods to analyze the regularization of 
a planar PWS fold bifurcation. 

In this paper, we regularize planar codimension one two-fold singularities that occur as the result of 
collisions of folds (quadratic tangencies) in both X~ and X + . We seek to identify PWS bifurcations as 
smooth bifurcations through regularization. We will study the fate of singular canards, pseudo-equilibria 
and limit cycles that can occur in our PWS system. We illustrate our analytical results with numerical 
simulations and show how the regularized system can undergo a canard explosion phenomenon. 

The paper is organized as follows. In section [2] we set up the problem, define the two-fold singularities 
we wish to regularize and present our PWS planar system in a normalized form such that the sliding regions 
retain their character under parameter variation. In section [3] we describe those properties of the PWS 
system that we wish to regularize, paying particular attention to singular canards, pseudo-equilibria and 
limit cycles. Then in section [4] we present a regularized version of our PWS system, using the approach 
of Sotomayor and Teixeira [301 . Before beginning our analysis, we collect together all our main results in 
section [5] giving the reader a concise summary of what is to come. In section [6] we carry out a blowup 
analysis m and show how singular canards persist the regularization. Our PWS system can also exhibit 
pseudo-equilibria, so in section [7] we consider how these unique PWS phenomena survive regularization. In 
section[Hl we show how limit cycles that are present in the original PWS system behave when regularized. In 
addition, we show how regularization can create another type of limit cycle that does not appear to be present 
in the original PWS system. For both types of limit cycle, we show how that the criticality of the Hopf 
bifurcation that gives rise to periodic orbits is strongly dependent on the precise form of the regularization. 
Some numerical results are presented in section [9] to illustrate our analysis. Our conclusions are presented 
in section [TO1 

2. Preliminaries. In this section we set up the problem, define two-fold singularities and present our 
PWS planar system in a suitable normalized form. Let x = (x,y) £ R 2 , /i £ R. Consider an open set 
(x, y) £ U x I and a smooth function / M (x) having 0 as a regular value for all y £ X. Then E C U defined 
by E = / / 7 1 (0) is a smooth ID manifold. The manifold E is our switching boundary. It separates the set 
E + = {( x,y ) £ Ulf^x.y) > 0} from the set E_ = {(x,y)|/ M (x, y) < 0}. We introduce local coordinates so 
that /^(x, y) = y and so E = {(x, y) £ U\y = 0}. From now on, we suppress the subscript y. 

We consider two smooth vector-fields X + and X~ that are smooth on E + and E_, respectively, and 
define the PWS vector-field X = (X~,X + ) by 


X(x,/T) 


X (x,/it) for x € E_ 
X + (x,y) for x £ E + 


( 2 . 1 ) 


Then 

• E cr C E is the crossing region where (A' + /(x, 0, y)(X~ f(x, 0, y)) = X^(x, 0, y)X%(x, 0, y) > 0. 


1 For example, they considered crossing, stable and unstable sliding, pseudo-saddle-nodes and two-folds. 
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• T, s i C E is the sliding region where (X + f(x, 0, p)){X~ f{x, 0, /z)) = X^{x, 0, y)Xf (x, 0, /i) < 0. 
Here X ± f(-,p) = V/ • X ± (-,p) denotes the Lie-derivative of / along X ± (-,p). Since f(x,y) = y in our 
coordinates we have that X±f = X±. 

In the sliding region, the vector fields on either side of E s i point either toward or away from E s /. In this 
case, in order to have a solution to our system in forward or backward time, we need to define a vector field 
on E s i. There are many possibilities, depending on the problem being considered. One of the most widely 
adopted definitions is the Filippov convention m, in which the sliding vector field X s i(x,fi) is taken to be 
the convex combination of X + and X ~: 

X s i(x.,fj,) = aX + (pc,n) + (1 - cr)X“(x, n), (2.2) 

where a £ (0,1) is such that X s i(x,p) is tangent to E s ;. In this case, 

= X ~f(x, 0,n) 

X-f{x,0,ri-X+f(x,0,liY 

The sliding vector field X s ;(x,/x) can have equilibria ( pseudo-equilibria , or sometimes quasi-equilibria 
HU). Unlike in smooth systems, it is possible for trajectories to reach these pseudo-equilibria in finite time. 
An orbit of a PWS system can be made up of a concatenation of arcs from E and E±. 

2.1. Two-fold singularities. The boundaries of E s ; and E cr where X + f = X£ = 0 or X~ f = Xf = 0 
are singularities called tangencies. The simplest tangency is the fold singularity, which is defined as follows. 

Definition 2.1. A point q £ E for p £ X is a fold singularity if 


X + f(q,»)= 0, X+(X+f)(q,y,)^ 0, 

(2.4) 

or if 


X~f{q,n) = 0, X-(X~f)(q,ia)^ 0. 

(2.5) 

A fold singularity q with X ± f(q,p) = 0 zs visible if 


X ± (X ± f)(q,y)^ 0, 

(2.6) 

and invisible if 


X ± (X ± f)(q,p)^ 0. 

(2.7) 


Note that, for y sufficiently small, the inequalities in (12.41) and (12.51) are equivalent to the following 

d x X+{q,0)X+(q,0) ± 0, d x X^(q,0)Xf (q,0) ^ 0. (2.8) 

In this paper, we consider the case of the two-fold singularity, when there is a fold singularity in both 
X±. In particular, we suppose that X± have tangencies at q ± = q^ip) £ E, respectively, which collide for 
p, = 0 at q = Q ± (0) with non-zero velocity. Hence ( q + — q~)'{ 0) ^ 0. 

Definition 2.2. We say that the two-fold singularity q is 

• visible if q + and q~ are both visible; 

• visible-invisible if q + (q~) is visible and q~ (q + ) is invisible; 

• invisible if q + and q~ are both invisible. 

The three different types of two-fold singularity are shown in Fig. 12.11 In the case of a single fold 
singularity, it is known that both the visible and invisible cases are structurally stable [T5] p. 232]. The 
regularization of the visible case was studied in [29] • Filippov m Figs. 58, 59] also considered the case 
of a single cusp singularity, which can be either visible or invisible. The cusp singularity is known to be 
structurally unstable, bifurcating into two tangencies PH Figs. 76, 77], which are on the same side of E. 
Kuznetsov et al. [22] Fig. 9] considered these bifurcations, which they label DT 12 , together with the cases 
we consider here. But we feel that the cusp singularity is best left for future work, as part of the wider picture 
that includes cusp-fold and two-cusp singularities. The two-fold singularities that we consider are shown in 
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Figure 2.1. The three different types of two-fold singularity studied in this paper: (from L to R) visible, visible-invisible, 
invisible. Following Filippov m, we show neither flow directions nor any sliding vector field. 


Filippov m Figs. 64, 65, 67, 68], where they are termed type 3 singularities There are 7 different generic 
cases. These were subsequently called VVlj, V/i_ 3 ,//i ,2 by Kuznetsov et al. [22]; a notation that we will 
find useful to adopt. Other authors mm refer to two-fold singularities as fold-fold singularities , which can 
be hyperbolic (visible), elliptic(al) (invisible) or parabolic (visible-invisible). Two-folds in R 3 were considered 
by the present authors in m- 

2.2. Normalized equations. In this section, we derive a normalized form for the equations near a 
two-fold singularity at (q,p) = 0 in R 2 . By Taylor-expanding X±, we have, for y > 0, 

x = X+(0) + 0(x + y + p), 

y = d v x£{0)y + d x x£{0)x + 0{\{x,y)\ 2 + p(x + y)), 

and, for y < 0, 

x = X-(0) + 0(x + y + p), 

y = d v x 2 (°)y + d * x 2 (°)(a ; - t l ) + 0(\(x - p,y) I 2 + p(x- p + y)). 

We now introduce x and t where 


x = 


t = 


X+(0) 


dyX+(0) 


sign(d y X+(0)) 


-t, 


^|X+(0)^X+(0)| 

which is well-defined, by virtue of (12.81) . Then, on dropping tildes, we have, for y > 0: 

x = 6 + 0(x + y + p), 
y = x + <D(y + px + x 2 ), 

and, for y < 0, 

x = a + 0{x + y + p), 

V = -P{x - p) + 0(y + p{x - p) + (x - p) 2 ), 
where <5 = sign(X^(0)d x X^"(0)) = ±1. The constants 


a = sign(9 y X^ (0)) 
13 = sign(d ?/ X 2 + (0)) 


AY(0) 

|X+(0)|’ 

+ ®y X 2 (0) 


\X+(0)d x X+ 


(2.9) 


( 2 . 10 ) 


2 Other type 3 singularities, shown in ESI Figs. 66, 69, 70, 71], have codimension greater than one (see |151 p. 239]). They 
include cusp-fold and two-cusp singularities. 
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are non-zero by (12.8|) . Later on, we will need to include higher order terms in our analysis. We introduce 
the following coefficients: 


C ± , X ± , and r] ± , 

so that (EH) becomes for y > 0: 

x = 5 + £ +a; + ^{x 2 + y + y), 
y = x + r/ + x 2 + x + V + 0(xy + y(x + y) + a; 3 ), 


( 2 . 11 ) 


( 2 . 12 ) 


and (12.101) becomes for y < 0: 

x = a + (~x + 0( x 2 + y + y), (2-13) 

V = - y) + ri~(x - y) 2 + x~V + 0(xy + yy + y(x - y) + (x - y) 3 ). 

Remark 2.3. De Carvalho and Tonon have given normal forms for codimension one planar P WS 
vector field^f |. However, we need ( f 2.9)) . h2.10\) and \2.12\) . \2.1S\) in this form in order to unfold several of 
the phenomena studied in this paper. 

The sliding vector field (12.21) is given by 

x = crXf (x, 0, y) + (1 — cr)Xf (a:, 0, y), (2.14) 

V = 0 , 

where a, defined in (12.31) . is given by 

_(-/3 + Q(x + y)){x- y) _ 

(—0 + O(x + y))(x — y)—(l + O(x + y))x 

The denominator in (12.151) is positive for stable sliding and negative for unstable sliding Sji. So if we 
multiply (12.141) by the modulus of this denominator, \(—j3 + 0(x+y))(x — y) — (l + 0(x + y))x\, corresponding 
to a transformation of time, we find on y = 0 that, in Sj, 

x = ±(—/3 + 0(x + y))(x — y)(5 + 0(x + y)) ± (—1 + 0(x + y))x(a + 0(x + y)), (2.16) 

y = 0, 


Equilibria of (12.161) are pseudo-equilibria, which we will study in section |TT~2l below. 

Within for y = 0 we find from (12.161) that 

x = {—/3S + 0(x))x — (a + 0(x))x = —{/36 + a)x + 0{x 2 ). 

Proposition 2.4. The fold q + = (0,0) is visible (invisible) from above if S = 1 (5 = —1), whereas the 
fold q~ = (y,0) is visible (invisible) from below if a/3 > 0 (a/3 < Q). Hence the two-fold q = (0,0) for y = 0 
is 

• visible if 8 = 1 and a/3 > 0; 

• visible-invisible if 5 = 1 (6 = —1) and a/3 < 0 (a/3 > 0); 

• invisible if 5 = — 1 and a(3 < 0. 

We also have that 


E s; : (3x(x -y)> 0, (2.17) 

E cr : f3x(x — y) < 0, (2-18) 


for x and y sufficiently small. The subset E“ z = E s ; D {aa < 0} of E s ; is the stable sliding region whereas the 
subset E^ = E s ; D {x > 0} of E s ; is the unstable sliding region. The subset T,~ r = E cr fl {cc < 0} of E cr is 
crossing downwards whereas the subset E+ = E cr . fl {x > 0} of E cr is crossing upwards. 

3 De Carvalho and Tonon (private communication) have indicated that they intend to publish a corrigendum to this paper, 
since their normal forms, as currently stated [7], can not distinguish between all the different planar PWS singularities. 
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Proof. These statements follow from simple computations. For example, to obtain the last part, we note 

that 


X£(x, 0, p)X 2 (x, 0, /jl) = (1 + 0(x + n))x(—/3 + 0(x + p)){x — /x), 

and use the definition of £ s j in section [2] together with (12.81) . □ 

Henceforth, in the visible-invisible case, without loss of generality, we will focus on the case 5 = 1, a/3 < 0 
so that the fold q + = (0, 0) is visible from above and q~ = (p, 0) is invisible from below (as in Fig. 12.11) . 

Since we perform a local analysis, we restrict attention to x and p sufficiently small so that statements 
(12.171) and (12.18[) in Proposition 12.41 about £ s ; and E cr , apply. The advantage of the form (12.91) and (12.101) 
of the normalized equations is that the sliding regions retain their character (stable or unstable) under 
parameter variation. 

For later convenience we introduce the parameter 52 defined by 

52 = /38 + a. 

To conclude this section, we state the following assumptions, which we make throughout the rest of the 
paper. 

Assumption 1. <5^0. 

Assumption 2. a/3 ^ 0. 

Assumption 3. 52 ^ 0. 

Assumption [l] and Assumption [2] are the normalized form of (12. 81) . when combined with (12.91) and (12.101) . 
In fact we have already set 5 = ±1. The significance of Assumption [3] will be explained in Remark 13. II below. 

3. Analysis of the PWS system. In this section, we analyze the planar PWS system (12.91) and (12.101) . 
together with (12.141) and (|2.15l) whenever we have sliding. We pay particular attention to singular canards, 
pseudo-equilibria and limit cycles that can occur in our system. The fold at q + = (0, 0) is fixed, whereas 
the fold at q~ = (p, 0) varies with p such that the two-fold at p = 0 bifurcates. Both pseudo-equilibria and 
sliding sections can appear, disappear or change character depending on whether q ± are visible or invisible. 
In addition, some of the two-folds at p. = 0 possess singular canards, which disappear for p ^ 0, and at least 
one two-fold can have a limit cycle. 

3.1. Singular canards. Trajectories can go from the attracting sliding region EJ ; to the repelling 
sliding region Ej), or vice versa, for p = 0. These trajectories, which we call singular canards [15], resemble 
canards in slow-fast systems [I]. A singular canard is called a vrai singular canard if it goes from the 
attracting sliding region E~ ; to the repelling sliding region in forward time. Singular canards that go 
from the repelling sliding region EJ, to the attracting sliding region £“ ; are called faux singular canards. 
Singular canards can only exist for p = 0 and when there is sliding in both x < 0 and x > 0. From (|2.17[) . 
we see that singular canards can only exist for /3 > 0. 

For the existence of singular canards in our PWS system, it is important to note that, in terms of the 
original time used in (12.141) . the two-fold on E s i can be reached in finite time. A simple calculation using 
L’Hopital’s rule shows that on E a j, for p = 0, 


lim x = (1 + /3) 1 fi. (3.1) 

x —^0 

There is no singularity at 1 + /3 = 0 since we need ff > 0 for sliding. So, by Assumption [3] we have a 
finite value of x on E s i for x —> 0 when p = 0. Hence it is possible to pass in finite time through x = 0 (the 
point separating attracting and repelling sliding regions, if they exist) at a two-fold. 

Remark 3.1. The case Q = 0 is degenerate, since lim x _,.o x vanishes. Geometrically this case corre¬ 
sponds to the linearized trajectories of X± having the same gradient on E s j. We shall not consider this case 
further (cf. Assumption^ above). 

Hence by m we conclude that singular canards exist in our PWS system. To decide whether they are 
vrai singular canards or faux singular canards, we need to consider the sign of x in (EH). We collect the 
results in the following proposition: 

Proposition 3.2. Singular canards in our PWS system exist if and only if f3 > 0. If fl > 0 (Cl < 0) 
then the singular canard is a vrai (faux) singular canard. 
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pseudo-equilibrium 


(d 7 ^ 0) 


no pseudo-equilibrium 


(v ~f~ 0 ) 



sliding 
(l* = °) 


crossing 

(d = o) 


Figure 3.1. The different types of two-fold singularities. Sliding is indicated by triple headed arrows. Cases //i ? 2 occur 
for 6 = — 1, the other cases for 5=1. 


One of the main objectives of this work is to establish persistence results of these singular canards under 
regularization. We will focus primarily on the persistence of vrai singular canards. 

The different types of two-fold, together with their flow directions and any sliding regions are shown 
in Fig. 13.11 Note that the visible two-folds W 12 and the visible-invisible two-folds VI 1-3 exist for 5 = 1, 
whereas the invisible two-folds 1 / 1,2 exist for 6 = —1. 

3.2. Pseudo-equilibria. As mentioned in section [3 the sliding vector field X s ;(x, /z) itself can have 
equilibria, called pseudo- equilibri^. These pseudo-equilibria are not necessarily equilibria of X ± (x, /z). In¬ 
stead they correspond to the case when X ± (x, /z) on S s i are linearly dependent. Filippov ([15] , P- 218) terms 
them type 1 singularities. They comprise three distinct topological classes; a pseudo-node, a pseudo-saddle 
and a pseudo-saddle-node. The following proposition describes the existence of pseudo-equilibria in (12.161) . 

Proposition 3.3. If 


a6 < 0 , 


(3.2) 


4 Hence there can be no pseudo-equilibria without a sliding vector field. 


























(3-3) 


then, for p sufficiently small, there exists a pseudo-equilibrium of \2.16\) at (x,y) = ( x ps ,0), where 

x ps = 05fl~ 1 p + ^(m 2 )- 
Also if 05fl~ 1 p < 0 then (x ps ,0) £ Ej ; and 

• for £2 < 0: (a: ps ,0) is a pseudo-saddle with local repelling manifold coinciding with EJ ; ; 

• for £1 > 0: (a; ps ,0) is a attracting pseudo-node. 

If 05fl~ 1 p > 0 then (x ps ,0) £ Ej] and 

• for £1 < 0: (£ ps ,0) is a pseudo-saddle with local attracting manifold coinciding with E^; 

• for £1 > 0: (:r ps ,0) is a repelling pseudo-node. 

If aS > 0, then ( x ps ,0) is not a pseudo-equilibrium. 

Remark 3.4. From Assumption^ we do not consider pseudo-saddle-nodes in our system. 
Proof. To find pseudo-equilibria, we set x = 0 in (12.1611 to get 

{—08 + 0{p + x))x + {—a + 0{p + x))x — {—05 + 0{p))p = 0, 


or 


(1 + 0 1 8a + 0{p + x))x = (1 + 0{p))p. 

Here we have used that 6 = ±1. Note that 

1 + 0~ 1 5a = 0~ 1 8~ 1 Pl ^ 0, 

by assumption. We can therefore solve this equation by the implicit function theorem to obtain 

x = x ps = 05 £2 -1 /i + 0{p 2 ). 

This is a pseudo-equilibrium if and only if (x ps ,0) £ E s j. We determine E s i as follows. Consider first 0 > 0. 
Then we have 


x < 0 

or 

X > p 

for 

M > 0, 

x = 0 

for 

P = 0, 



x < p 

or 

x > 0 

for 

p < 0 . 


Thus E s ; : sign(/z)x ^ (0, |/z|). Then since 

sign {n)x ps = 05n~ 1 \p\ + O{p0), 

we conclude that (ir ps ,0) £ E s ; for p sufficiently small, provided 05PL~ 1 = (1 + 0~ l 5a )~ 1 ^ (0,1). Since 
0 > 0, this condition is equivalent to a5 < 0. If, on the other hand 05Pl~ l = (1 + 0~ 1 5a)~ 1 £ (0,1), then 
a5 > 0, for 0 > 0 and so (x ps , 0) is not a pseudo-equilibrium. 

Next consider 0 < 0. Then 


E s ; : sign{p)x £ (0, \p\). 

Hence (a; ps ,0) £ E s ; for p 0 sufficiently small, provided 0 <5£I -1 = (1 + 0~ l 5a)~ 1 £ (0,1). For 0 < 0, this 
is equivalent to aS < 0. If, on the other hand 05Pl~ l = (1 + 0~ 1 5a)~ 1 ^ (0,1), then aS > 0, for 0 < 0 and 
so {x ps , 0) is not a pseudo-equilibrium. 

We conclude that (cc ps ,0) is a pseudo-equilibrium, where x ps is defined in (13.31) . provided (13.21) holds. 

If Xp S < 0 then ( x ps , 0) £ E , the region of stable sliding, and if x ps > 0 then (x ps , 0) £ E^, the region 
of unstable sliding (see Proposition ^. 41) . Combining this with the linearization of (12.161) about (cc ps ,0) gives 
the statements about stability. □ 

These results are summarized in Fig. 13.21 As mentioned earlier, two-fold singularities occur in (EH) 
and (12.101) when p, = 0. From Proposition 13.31 it follows that the two-fold singularity can be accompanied 
by significant changes in the nature of pseudo-equilibria around p = 0. For example, if £2 < 0 then x ps is 
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j36£l 1 /x < 0 


n < 0 

n > 0 


I 

Hi 

H 

^ t 

Co 

X 

= 0 

n < 0 


n > 0 


Xps 

X Xp S 

X 

psn V > 0 



n < 0 

Xps 

n>o 

Xps 




r 


x 


Figure 3.2. Results of Provosition UT3\ Here aS < 0 and x ps = fiSQ O(pP). 

a pseudo-saddle for all /j, ^ 0. The difference between /j, > 0 and /r < 0 is, where /3<5H _1 /r < 0, the pseudo¬ 
saddle is in which coincides with the associated repelling manifold of x ps . But, where /3<5fl -1 /r > 0, the 
pseudo-saddle is in Ejj, which coincides with the associated attracting manifold of x ps . Hence the attracting 
and repelling directions “switch” on passage through the two-fold at /r = 0. We can see this behaviour in 
[221 Fig. 10], where E s ; is the attracting manifold of the pseudo-saddle of VV 2 on one side of the two-fold, 
whereas it is the repelling manifold on the other. Similarly, when fl > 0 a pseudo-equilibrium goes from 
being attracting on one side to repelling on the other side. We see this behaviour for example in [22] Fig. 
11], for the VI 3 two-fold, where a repelling pseudo-node becomes an attracting pseudo-node. 

Another main objective of this paper is to understand how the behaviour of these pseudo-equilibria is 
modified when our governing equations are regularized. 

3.3. Limit cycles. A further phenomenon in the two-fold singularity is the existence of limit cycles. 
For our PWS system, it is clear from Fig. 13.11 (see also [22J Fig. 12]) that (local) limit cycles can occur in 
the invisible case II 2 where: 

5 = —1, a > 0, /3 < 0. (3.4) 

To study periodic orbits in this case, one can introduce a Poincare map Pq which takes {(x, 0)|x > 0} into 
{(x, 0)|x > /j,} under the forward flow of A' + and A' - (see Fig. l3.3l and [T5I p. 236]). The map P 0 is composed 
of {Tq" : xq 1 ^ Xi, the mapping from (xq,0), xq > 0, to (xi(xq),0), Xi(xq) < 0, under the forward flow of 
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Figure 3.3. The Poincare map Pq = a Q o ctq associated with the invisible two-fold singularity II 2 , where 5 = — 1, a > 0 
and (3 < 0. 


X + , and < 7 q : x\ K>• £ 2 ) the mapping from (xi,0), Xi < /x, to (* 2 (^ 1 )) 0), a^O^i) > Ab under the forward flow 
of X - . Clearly Po is only defined for those Xo for which X\ < fi. For X~ to map x\ < /.1 into X 2 > p we 
need x = a + C7 (.t + y + p) >0. Therefore a > 0. Hence we have the following lemma. 

Lemma 3.5. Consider jX^p . Then 

a o ( x o ) = -xo + A + x q + 0(xo(n + * 0 )), 

= —xi + 2 p + A~ x\ + (D(xi (/x + if)), 

where 

A ~ = -J/j ( aT?_ + + ( 3 - 5 ) 

H+ = -^( ?? ++C + +x + )- 

Proof. See [T5] p. 236]. □ 

The Poincare mapping Po = (Tq 0 <Tq therefore takes the following form: 

Po{x 0 ) = ctq (<t+ ( xq )) = xo + 2/x + ( A~ - A + )xl + Q(xl + pxo), (3.6) 

for those xo which satisfy the inequality: 

(To (xo) = -xo + A + x 0 + 0(x 0 + xop) < /x. (3.7) 


Proposition 3.6. Consider &3.4\ ) and suppose that 

A Il2 =A~~A+ = ^ (a( V ~ + Pv + ) + 0(C + x ~ + a(C + + X + ))) , (3-8) 

is non-zero. Then, for /xA Jj < 0 sufficiently small, the PWS system has a family of periodic orbits. These 
periodic orbits correspond to fixed points of Pq of the following form 

xo ( n ) = y /— 2 fj , Ajj 3 + O ( p ). (3.9) 

The periodic orbits are attracting for A jj 2 < 0 and repelling for A jj 2 > 0. 

Proof. We obtain (13.91) by solving the equation P 0 (;ro) = xo for Xq > 0 using the implicit function 
theorem. For /x small, Xo (/x) is positive and satisfies the inequality (13.71) . For stability we compute the 
derivative of (13.61) to get: 


Pq(xq (p)) = 1 + 2A// 2 x 0 (/x) + O(p). 
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Figure 3.4. In (a): the VI 3 PWS two-fold bifurcation and in (b): the classical curved critical manifold e.g. appearing in 
the slow-fast van der Pol system. In (b) S a and S r denote the attracting and repelling part of the critical manifold. Note also 
the different use of triple-headed arrows. In (a) triple-headed arrows are used to indicate sliding whereas they in (b) are used 
to indicate fast orbit segments described by a set of layer equations. The thick curves illustrate singular closed curves which 
for case (a) will be become periodic orbits upon regularization. 


If A// 2 < 0 (A u 2 > 0) then 0 < P'(0) < 1 (P'( 0) > 1) for p, sufficiently small. The result follows. 0 

Finally in this subsection, we note that the visible-invisible case VI3 does not appear to have limit 
cycles. However, by straightening out the flow within E _ , so that E becomes curved and quadratic at the 
fold, the VI3 case clearly resembles the classical slow-fast curved critical manifold with singular cycles. This 
similarity can be seen in Fig. 13.41 where we illustrate (a) the VI3 two-fold and (b) the slow-fast equivalent 
as it appears, for example, in the van der Pol system (see also [21] Fig. 5]). There are two singular canard 
cycles in Fig. 13.41 (b), shown as thick curves, each composed of fast segments (with triple-headed arrows) 
and slow segments on the curved slow manifold. For e sufficiently small, this slow-fast system is known 
to produce a canard explosion phenomenon [21j . in which the singular cycles become limits of a family of 
rapidly increasing periodic orbits as e —> 0. However, when considering Fig. 13.41 it is important to highlight 
that there is no time scale separation in the PWS system. Nevertheless, we shall see that the regularized 
PWS system possesses a hidden slow-fast structure near the discontinuity set. This will allow us identify 
the singular cycles in Fig. 13.41 (a) as limits of periodic orbits of the regularization, and hence strengthen 
the connection between (a) and (b) further. The singular cycles in Fig. 13.41 (a) are also illustrated as thick 
curves, but in comparison to the cycles in (b), they are composed of an orbit segment of X ~, within E _ , 
and a sliding segment on E. 

Another objective of this paper to understand the existence of limit cycles under regularization. Further 
details are presented in section [8] 

3.4. Summary of two-fold properties. We conclude this section with Table [3T7TI which summarizes 
properties of the seven two-folds, for future reference. 


Table 3.1: 1: type of two-fold singularity [visible (VV), visible- 
invisible (VI) and invisible (II), following [22]]. 2: value of 5. 3: 
sign of a. 4: sign of /3. 5: sign of fi. 6: type of singular canard 
(x=no canard). 7: type of pseudo-equilibrium after bifurcation 
(PS=pseudo-saddle, PN=pseudo-node, x=no pseudo-equilibrium). 
8: possibility of limit cycle. 


1 

2 

3 

4 

5 

6 

7 

8 

Two-fold type 


a 

p 

n 

Singular canard 

Pseudo-equilibrium 

Limit cycle 

VVl 

1 

+ 

+ 

+ 

vrai 

X 

no 

VV 2 

1 

- 

- 

- 

X 

PS 

no 

VI\ 

1 

+ 

- 

± 

X 

X 

no 

VI 2 

1 

- 

+ 

- 

faux 

PS 

no 
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Vh 

1 

- 

+ 

+ 

vrai 

PN 

jj 

III 

-1 

- 

+ 

- 

faux 

X 

no 

Ih 

-1 

+ 

- 

+ 

X 

PN 

yes 


4. Regularization. It is natural to ask how the results in section [3] are affected by regularization. 
For example, there is the question of the persistence of the singular canards isection ld.il) . Then, of course, 
pseudo-equilibria (section 13. 21) can not exist in a smooth system, so we will need to look at the existence and 
behaviour of equilibria in the regularization. Finally there is the need to understand both the fate of the 
limit cycles (section 13.31) in the invisible two-fold 11-2 under regularization and the possibility of other limit 
cycles which may appear in the regularized system. 

There is a number of ways that the original PWS vector field A' = (X“,AT + ) can be regularized. We 
follow the approach of Sotomayor and Teixeira [32]. We define a C fc -function (1 < k < oo) cf>(s) which 
satisfies: 


where 




1 for s > 1, 

S(-l,l) for se(-fl), 
— 1 for s < — 1, 


(4.1) 


<j) (s) > 0 within s £ (—1,1). 


(4.2) 


Then for e > 0, the regularized vector-field X £ (x, y) is then given by 


^e(x,fi) = ^I + (x,/i)(l + ^(e 1 y)) + ix (x, y)(l — (j)(e l y)). 

Note that 


(4.3) 


-Xe(x, y) = A ± (x, y) for y^±e. (4.4) 

The region y £ (—e, e) is the region of regularization. The system (14.31) has a hidden slow-fast stucture which 
is unfolded by the re-scaling (see also [23] 1 

V = e”V (4.5) 

In terms of this re-scaling the region of regularization becomes y £ (—1,1). Then, using (12.91) and (12.101) . 
the regularized system (POl) becomes 

x = \ ((<5 + 0(x + y + ey))( 1 + + (a + 0{x + y + ey))( 1 - <j>{y ))), (4.6) 

ey = \{(x + <D(yx + x 2 + ey))( 1 + (j)(y)) + (~0(x - y) + 0(y{x - y) + (x - yf + ey)){ 1 - </>(£))) . 

This is the slow system of a slow-fast system. The y variable is fast with (D(e -1 ) velocities and x is the slow 
variable with 0(1) velocities. The limit e = 0: 

±= \ i( S + °( x + m))(4 + 4>{y)) + (a + 0{ x + n))(l - 4>{y ))), (4.7) 

o = \ {(x + 0{yx + x 2 ))(l + (j){y)) + (-/3(x - y) + 0(y(x - n) + (x - /x) 2 ))(l - <t>(y))) ■ 

is called the reduced problem. When we re-scale time according to r = -^-t, (POll becomes: 

x' = e ((5 + 0(x + p + ey))( 1 + <j>(y)) + (a + 0(x + p + ey))( 1 - <j>{y ))), (4.8) 

y’ = (x + 0(yx + x 2 + ey))( 1 + 4>{;y)) + (-/3(x - y) + 0(y(x - y) + (x - y) 2 + ey)){ 1 - <f(y)), 


5 The regularized two-fold VI 3 does possess a limit cycle. See Theorem 18.81 
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with ()' = -4-. This is the fast system of a slow-fast system. The limit e = 0: 

x' = 0, (4.9) 

y' = (x + 0(px + a; 2 ))(l + <j>(y)) + {-/3{x - p) + 0(p(x - p) + (x - m) 2 ))0 - </>($))■ 


is called the layer problem. Time r is the fast time and time t is the slow time. 

Remark 4.1. The re-scaling used in (O is necessary to identify the slow-fast structure hidden in 
O- It is not due to loss of hyperbolicity and so the use of ra is different from the scaling used later in 
section^ in connection with the blowup of a nonhyperbolic line. We use m to cover y £ (—1,1). Outside 
this region we will use o and the exact PWS analysis to describe the dynamics. 

In J23i \24\ \25f the re-scaling is replaced by a transformation y = ry , e = re, (y,e) £ S 1 , making 

it possible to illustrate, in one diagram , the dynamics of outside the region of regularization and the 
slow-fast dynamics within (by inserting the cylinder ( x , ( y,e )) glxS 1 at y = 0). However it does nothing 
to address any loss of hyperbolicity that may occur. 

In this paper, we apply and extend Fenichel theory m m m nu of singular perturbations to study 
these regularized equations (I4.6D - (I4.9D . allowing us to go from a description of the singular limit e = 0 to a 
description for e > 0. This approach has the advantage that it is geometric so, for example, we are able to 
solve persistence problems by invoking transversality. It will also be important to identify the singular limit 
e = 0 with the original PWS system. 

The key to the subsequent analysis in this paper is the following result (a related result is given in [24[ 
Theorem 1.1], in a slightly different form): 

Theorem 4.2. There exist critical manifolds S ajT and q, given by 


, 1 - <j)(y) _ x + 0(x 2 + fix) 

1 + <t>(y) /3(x-p) + (D{p(x-p) + (x~ p) 2 )’ 

, 1 - 4>(y) _ x + 0(x 2 + px) 

l + (f{y) P{x - p) + 0(jj,(x - p) + (a: - p) 2 ) ’ 
q: x = 0 , p = 0 , y £ (- 1 , 1 ), 


(4.10) 


(4.11) 


which are normally attracting, normally repelling and normally non-hyperbolic, respectively. On the critical 
manifolds ,S a ^ r , the motion of the slow variable x is described by the reduced problem 0, which coincides 
with the sliding equations \2.1Iff and \2.15\) . 

Proof. This follows from simple computations. Recall that the critical manifolds are fixed points of the 
layer equations (Oil . Each manifold is normally hyperbolic if the linearization has one non-zero eigenvalue. 
Otherwise it is non-hyperbolic. 0 

Remark 4.3. The non-hyperbolic line q of the regularized system corresponds to the two-fold q of the 
PWS system. 

As in our previous paper |T8] it is useful to introduce the following function w = w(y ): 


i - <Kv) 




i + 4>(y) 


(4.12) 


which appears on the left hand side of (14.101) . since this reduces the complexity of subsequent expressions. 
For y £ (—1,1) we have w £ (0, oo). Also 


w'{y) = 


Wiy) 

Y+m? 


< 0 , 


within y £ (—1,1), and the critical manifolds S air are therefore graphs of y over 


S r ■ 


y = w 


y = w 


x + 0(x 2 + px) 
ft(x — p) + <D(p(x — p) + (x 
x + 0 (x 2 + px) 
/3(x — p) + 0(p(x — p) + (x 



(z,0) e E s ,, 
0,0) e E+, 


(4.13) 
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where 

w~ 1 {z) = cj)- 1 (-j-Tpr) ’ ^t 0 ’ 00 ) 

is the inverse of w. 

In the sequel, we will also need the extended problem , which is given by 

x' = e ((<5 + C + £ + 0(x 2 + p + ey))( 1 + < p(y)) + {a + (~x + 0(x 2 + p + ey))( 1 - <p(y))) , (4.14) 

y' = (x + r/ + x 2 + ex~y + 0 (exy + p(x + ey) + z 3 ))(l + <j>(y)) 

+ (~P(x -n)+rf{x- p ) 2 + e\~y + 0 (exy + pey + p(x - p) + {x - y) 3 )){l - 0 {y)), 
e' = 0, 

p! = 0 . 

Here we have used (12.121) and (12.131) . 

Finally, when studying canards in section [6l it will also be useful to introduce a new time f, defined by 

df = (1 + <j)(y))dT, (4.15) 

so that our extended problem (14.141) . now only defined within y £ (—1,1), becomes: 

x' = e [5 + C, + x + 0(x 2 + p + ey) + (a + (~x + 0(x 2 + p + ey))w(y)) , (4-16) 

y' = x + rj + x 2 + ex~y + 0(exy + p(x + ey) + x 3 ) 

+ {~P{x - p) + y~{x - p) 2 + X ~y + 0 (exy + pey + p{x - p) + (x - p) 3 ))w(y ), 
e' = 0, 

p = 0. 


where, by abuse of notation, we take ()' = in section \E\ only. Note that w(y) appears on the right hand 
sides of equations (14.161) . 

Remark 4.4. If we return to our original y variable using then the critical manifolds S a , r become 
graphs y = eh(x,p) which, for e = 0, collapse to S fl j. 

Remark 4.5. The critical manifolds S a , r and q are shown in Fig. \J h 1 \ for j3 > 0. It is clear that a 
singular canard exists for p = 0 only. 

4.1. Fenichel theory. Consider compact sets / ± contained within Ej), respectively. Then according 
to Fenichel theory mmmm, there exist slow manifolds of (14.31) of the form 




y = w 

y = w 


x + 0 (x 2 + px) 
/3(x — p) + 0(p(x — p) + (x 
x + 0 (x 2 + px) 
P(x — p) + 0 (p(x — p) + (x 



O(e), Or, 0) ere Ejj, 
O(e), (x,0)£l+cZ+. 


(4.17) 

(4.18) 


The flow on S 0tF and S r>e is therefore O(e)-close to the flow of the sliding equations. The potential non¬ 
uniqueness of S a ,e and S r , e only manifests itself in 0(e~ c/ ' e ) small deviations. Henceforth we shall fix copies 
of S a , e and S r . e . 

4.2. Loss of hyperbolicity. From Theorem l4.21 it follows that the analysis in section[3]can be directly 
applied to singular canards on the critical manifold for the limiting regularized system. In contrast, a 
(maximal) canard of the regularized system (14.81) appears as an intersection of the extension by the flow 
of the Fenichel slow manifolds S a ,e and SV,e to a vicinity of the non-hyperbolic line q. A canard of the 
regularized system is called a vrai (faux) canard if it goes from S aie (5 r , e ) to S r , e (S a , e ) in forward time. 

Remark 4.6. Our focus is mainly on vrai canards, usually referred to as canards \ 19 i 1 31 ^ . a convention 
we shall adopt. 

To study the extensions of S a , e and S r<e near q , we will use the blowup method of Dumortier and 
Roussarie ®mm, in the formulation of Krupa and Szmolyan |19] . 
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Figure 4.1. The critical manifolds S a ,r for different values of p together with the non-hyperbolic line q for f3 > 0. The 
critical manifolds agree with the sliding regions of the P WS smooth system. In particular, stable (unstable) sliding becomes an 
attracting (repelling) critical manifold of the regularization carrying slow dynamics described by the sliding equations. 


Consider the following sections: 


A T : x = TP, (4.19) 

where p is a small positive constant. Then Fenichel’s slow manifolds S ai£ and S r ^ e intersect A^, respectively, 
in 


y = w ®(P x ) + O a>r {p + p, T e), (4.20) 

from (14.171) and (14.181) . This will become useful later when extending S (lif and S r<e . 

We are also interested in studying the regularization of pseudo-equilibria and limit cycles. The non- 
hyperbolicity of the line q defined in (|4.11|) will again complicate the analysis and we will also use the blowup 
method to obtain an accurate and complete description of the regularization of these phenomena. 

5. Main results. In this section, we anticipate our main regularization results, in a form convenient 
for the reader. The connection between sections of the paper devoted to PWS results and to regularization 
results is shown in Table EU 


Table 5.1 

Connection between phenomena in the PWS and regularized systems 


PWS phenomenon 

Regularized phenomenon 

Singular canards (section [3T|) 

Canards (section [6]) 

Pseudo-equilibria (section 13.21) 

Equilibria (section 0) 

Limit cycles (section 13.31) 

Limit cycles (section [8]) 


5.1. Canards. From section I5TT1 singular canards of the PWS system can only exist for p = 0. From 
Proposition 13.21 singular canards exist for 


/3 > 0 and > 0. (5.1) 

Canards of the regularized system EM are considered in section [Gj The main result of that section is that 
singular canards survive the regularization in the following sense: 
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THEOREM 16.41 Assuming m, the regularized system El has a maximal canard at 

fz = v / eM 2 ,c(\/e) = 0(e), 

where p 2 , c is given by 16. 22\) . which is 0(yfe)-close to the singular canard. 

Remark 5.1. Canards are non-unique because the slow manifolds are non-unique. But often in the 
literature, any solution following 5 r , e for an 0(l)-distance is called a canard. Hence the canard in Theorem 
which is obtained geometrically as the transverse intersection of fixed copies of 5 0 , e and 5 r , e , is referred 
to as the maximal canard. 

5.2. Equilibria. From section EOl pseudo-equilibria of the PWS system can only exist for 

ad < 0. (5.2) 

Equilibria of the regularized system (1X31) are considered in section [3 We have two main results. Our first 
result is Theorem l7.51 which shows that, for 12 < 0, the equilibria are saddles and, for 12 > 0, there is a Hopf 
bifurcation at p = 0(y/e), after which the equilibria become foci, and then nodes, as p varies. 

theorem 17.51 Assuming \5.2 1) . the regularized system El has a smooth and locally unique family of 
equilibria 


(x, y) = (x*,y*)(p, y/e) = ^ + p- 1§a + 0 (m 2 + e),°( e )j 

where (x* (p,0),y* (p,0)) = (x*(p,0),0) agrees with the family of pseudo-equilibria for the PWS system (see 
Provosition \3.3\ ) and, in particular, d li x*(p,0) ^ 0. For 

12 < 0: The family of equilibria consists of saddles and does not undergo any bifurcation. 

12 > 0: The family of equilibria undergoes a Hopf bifurcation at p = y/ep- 2 ,H{y/t) = 0(e), where P2,h is 
given by {7.Iff) . The first Lyapunov coefficient is given by a = a 2 \fe + 0(e) where 02 , given by 
\7. 15\ ), depends upon the regularization function <f>. 

Our second result for equilibria of the regularized system (1X31) . Theorem ES is perhaps surprising, in 
that it shows that the criticality of the Hopf bifurcation depends on the regularization function <f>. 

theorem 17.61 Suppose that 12 > 0 so that, from Theorem \7.5\ there exists a Hopf bifurcation. Then for 
e sufficiently small, provided 


5{C + X~) ~ a{t + + X + ) ^ 0, CLSD 

where , x ^ are the higher order coefficients in \2.11\) . the first Lyapunov coefficient a can be positive, zero 

or negative, depending on the regularization function </>. 

5.3. Limit cycles. From section ITTTTTl limit cycles of the PWS system can exist for the case II 2 , which 
occurs when 


6 = —1, a > 0, 0 < 0. (5.3) 

Limit cycles of the regularized system (14.31) are considered in section [8] where, following Theorem 17.51 above, 
we find limit cycles provided: 


12 > 0 and ad < 0. (5-4) 

This leads to the regularization of II 2 and of VI 3 , which we denote by I/| and EJJ, respectively. We have 
two main results for limit cycles of the regularized system. The first main result Theorem 18.81 shows how 
small amplitude periodic orbits due the Hopf bifurcation in Theorem 17.51 can be connected to 0(1) (with 
respect to e) amplitude periodic orbits. 
theorem 18.81 For e sufficiently small: 

IIThere exists a C k -smooth family of locally unique periodic orbits of the regularized system that 
is due to the Hopf bifurcation in Theorem \ 7. 5\ If a 2 < 0 (02 > 0) where 02 is the first Lyapunov 
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coefficient as defined in \1.15\ ), then the periodic orbits are attracting (repelling) near the Hopf 
bifurcation. If 

A// 2 = (■ a(rj~ + (3r] + ) + ffiC + X~ + a{ C + + x + ))), (EUD 

is negative (positive) then the periodic orbits for Aj^/x < —cffie, c > 0 sufficiently large, are attract¬ 
ing (repelling). The periodic orbits for p = 0(1) (with respect to e) are continuously O(e)-close to 
the PWS periodic orbits in Provosition \3.6\ 

V JJ: There exists a C k -smooth family of small periodic orbits of the regularized system 14-3\ ) that is due 
to the Hopf bifurcation in Theorem There also exists a C k -smooth family of periodic orbits 

that are 0 ( 1 ) (with respect to e) in amplitude and which undergo a canard explosion, where the 
amplitude changes by 0(1) within an exponentially small parameter regime around the canard value 
/i = \fepi,c(\[e) (see Theorem B- If a -2 <0 (a 2 > 0) where a 2 is the first Lyapunov coefficient as 
defined in |7. 15\ ), then the periodic orbits are attracting (repelling) near the Hopf bifurcation. If 

Avh = 3a/3 + + ^ + 1 ' ) ^~ + d 8 - 11 !) 

is negative (positive) then the (D(l)-periodic orbits are attracting (repelling). 

We conjecture on the connection of the two families of periodic orbits in VJJ: 

Conjecture [T] The two families of periodic orbits in VJJ belong to the same family of locally unique 
periodic orbits. 

The second main result for limit cycles of the regularized system (14.31) shows how an open set of regu¬ 
larization functions can induce at least one saddle-node bifurcation in the periodic orbits of Theorem 18.81 

theorem 18. 1 81 Suvvose j 7 .19( 1 and j5.4\ ). Then for e sufficiently small: 

JJJ: There exists an open set of regularization functions such that the periodic orbits in Theorem \ 8 . 8 \ 
case JJJ, undergo at least one saddle-node bifurcation. 

V JJ; Suppose, in addition, that Conjecture [7] holds. Then there exists an open set of regularization func¬ 
tions such that the periodic orbits in Theorem Iff.ffl case VJJ, undergo at least one saddle-node 
bifurcation. 

In section [9j we illustrate the results in Theorem 18.181 by applying two different regularization functions 
to two model systems of JJJ and VJJ. The regularization functions are such that each case, JJJ or VJJ, will 
have a saddle-node bifurcation for only one of the regularization functions. 

6 . On the existence of canards. In section 13.11 Proposition 13.21 we showed that our PWS system 
(EH and (12. 1011 together with the sliding vector field (12.141) and (12.151) contains singular canards for /3 > 0. 
In this section, we aim to discover the fate of these singular canards in the regularized system. To do this, 
we focus on dynamics in the region of regularization y £ (—e, e), or y £ (— 1 ,1), as described by equations 
(14.161) . which are written in terms of the new time t, defined in (14.151) . We work with t in this section only. 

As discussed in section l4~2l Fenichel theory breaks down on the non hyperbolic line q, defined in (14.111) . 
We use the blowup method [annum to deal with this line. We introduce the quasi-homogeneous blowup , 
given by 

(x,e,p) = (r ai x,r a 2 e,r a 3 p), 

where the number r is called the exceptional divisor. By this transformation, the line q is blown up to a 
cylinder 

S 2 x (-1,1) = {((x,e,p.),y)\x 2 + e 2 + fi 2 = 1, y £ (-1,1)}. 

When r = 0, the blown-up coordinates collapse to the non-hyperbolic line q. 

The weights (ai, 02 , 03 ) are chosen so that the vector field written as a function of the blowup coordinates 
has a power of the exceptional divisor as a common factor. By transforming time using this common 
factor, it is then possible to remove the exceptional divisor and so de-trivialize the vector-field on q. By 
substituting the quasi-homogeneous blowup into (14.161) and removing the exceptional divisor, it turns out 
that (ai, a 2 , 03 ) = (1, 2,1). So we have the following blowup of q: 

x = rx, e = r 2 e, p = rp, 
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Figure 6.1. The two-sphere S 2 and the charts K 1 — 3 . 


with r > 0, (x, e,/x) £ S 2 . Note that this blowup does not depend on y. The new phase space is therefore 

(( x,e,p),r,y) eS 2 xI + x (-1,1). 

To describe the dynamics on the blowup space we consider the following charts: 


chart 

«i : 

X = 

-1, X = -n, e = r\e 1 , p = r 3 pi 

( 6 . 1 ) 

chart 

k 2 : 

e = 

1 , x = r 2 x 2 , e = r\, p = r 2 p 2 , 

( 6 . 2 ) 

chart 

H3 ■ 

x = 

1 , x = r 3 , e = rle 3 , p = r 3 p 3 . 

(6.3) 


The chart k 2 is called the scaling chart or family rescaling chart mmm- The charts ki j3 are called phase 
directional charts. The point (x, e,/z) = (0,0,0) has been blown up into the planes r% = 0 : i = 1,2,3. The 
two-sphere S 2 and the charts ki _ 3 are shown in Fig. 16.11 We adopt the convention that the subscript n of 
each quantity is used when we are working in chart n n . 

The following coordinate change between charts K\ and k 2 will be important in what follows: 


« 2 i : (»T,y,ei,^i) '-t (x 2 ,y,H 2 ,r 2 ) = (—e x 1/2 , y, e x 1/2 Pi, r ly /e[), (6.4) 

«i 2 : {x 2 ,y,r 2 ,p 2 ) (r 1 ,y,e 1 ,p 1 ) = (-r 2 x 2 , y, xf 2 ,-xf 1 p 2 ), (6.5) 

defined for x 2 < 0 and e\ > 0 respectively. The change between charts n 2 and n 3 is given by: 

K 32 : (x 2 ,y,r 2 ,H 2 ) >->■ (r 3 ,y,e 3 ,p 3 ) = (r 2 x 2 , y, xf 2 ,xf 1 p 2 ), ( 6 . 6 ) 

«23 : (r3,y,e 3 ,P3) (x 2 , y, P 2 , r 2 ) = (e 3 2 , y, e 3 ~ 1/2 y 3 , r 3 s/e 3 ), (6-7) 

defined for x 2 > 0 and £3 > 0 respectively. 


We now describe the dynamics in each chart, beginning with the chart k 2 . Sections 16.21 and 16.31 describe 
the dynamics in charts ki j3 . Finally, section [6741 combines the results from charts n\- 3 to prove Theorem l6.4l 
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6 . 1 . Chart k 2 . The extended problem (14.161) written in chart n 2 becomes 

± 2=8 + aw(y) + r 2 (C, + x 2 + Cx 2 w(y)) + 0(r 2 (y 2 + r 2 )), (6.8) 

y = x 2 + r 2 (y + xl + x + y) , 

+ (~P(x 2 ~ Hi) + r 2 (y~(x 2 - M 2 f + X~y)) w{y) + 0(r 2 (y , 2 + r 2 ))), 
f 2 = 0 , 

M 2 = 0, 

where we have multiplied time by r 2 . For j3 > 0 there is an invariant line l 2 given by: 

h- x 2 € R, y = y c , r 2 = 0, n 2 = 0, (6.9) 


where 

y c = w -\r 1 ) = r 1 (|^ 4 ) e (-1,1), 

carrying the special solution: 

x 2 = /? _1 nt 2 , y = y c , r 2 = 0, y 2 = 0, 

where t 2 is the time used in (16.81) . 

Consider the following sections 

AJ = {(X 2 ,y, r 2 ,y 2 )| x 2 = +v~ 1 }, 


( 6 . 10 ) 


( 6 . 11 ) 


where v is small and positive. The line l 2 intersects Aj in 

h n AJ : x 2 = T^ _1 , y = tic , r 2 = 0, y 2 = 0. (6-12) 

The line l 2 and the sections AJ are shown in Fig. 16.21 

6.2. Chart K\ . We now describe the dynamics in chart n\ and relate them to the dynamics in chart 
k 2 . The extended problem (14.161) written in chart K\ becomes 

h = — rieiFi(n,y, (6.13) 

y = -i + P(i + m )w(y) + 0{n), 

ei = 2elF 1 (r 1 ,y,e 1 ,y 1 ), 

Mi = ^HxiFi{ri,y,ei,yi), 
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where we have divided the vector-field by rq and set F\{r\,y , ei, Mi) = 5 + aw(y ) + 0(r\). The (D(ri)-terms 
include the constants . 77 , from H2.11D but they will not play a role in this section and we therefore 

suppress them. The line r± = 0, y = ui -1 , ei = 0 is a line of fixed points, provided /?(1 + Mi) > 0, 

since w(y) € (0, oo). In particular, this line includes the point r\ = 0, y = y c , ei = 0, mi = 0 because we 
focus on P > 0. There exists an attracting center manifold M a p of this line, given by: 

M - i: » = («TTS)) +0(ri + ei) 

and, within r i = 0, a center manifold C a p, given by: 


C a . 


l : y = w 


1 


/3(1 + Mi) 


0(e i). 


Within C a , i, there exists an invariant line: 


la, l : ri = 0, y = y c , ei > 0, Mi = 0. 

Recall (16.101) . The center manifold C^i is overflowing and hence unique near /„ i if > 0 for > 0. From 
(16.131) . we therefore have 


ei = 2 ej(6 + aw(y)) = 1 VL + 0(e 1 ), 

with n = 0 and w(y) = f3~ l + 0(e i) near l a ,\. Hence since /? > 0, ei > 0 for ei > 0 if H > 0. For fl < 0, we 
have ei < 0 and so C a , i is non-unique near l a ^ in this case. 

The manifold M„ -\ has invariant foliations, which we denote by M a , i(e), with e = r\e\ = const. The 
sub-nranifold M a ,\{e) intersects the section r\ = p , corresponding to A - (14.191) . in 

y = + Q{p + m + e) = y c + 0 (p + p + e), 

which agrees with (14.201) since p = r±pi = pp ,Hence M a , i(e) is the extension of S a , e into chart n\ near the 
line q (ignoring exponentially small terms). 

In order to relate the dynamics in chart k\ to the dynamics in chart « 2 , we use the coordinate change 
« 2 i in (16.51) . The section A 2 defined in (16.111) then becomes 

A(“ : ei = i/ 2 . (6.14) 

The manifold C a ,\ intersects A(" in 

C a ,i D A(~ : ei = v 2 , n = 0, y = y c + <D{v 2 + pi) 

and hence, by the conservation of e, we conclude that the intersection of M a i(e) with A]" is O(v / e)-close to 
Cap. 

The line l a , i and the section A(" are shown in Fig. 16.31 

6.3. Chart K 3 . The extended problem (14.161) written in chart K 3 becomes 


^3 = r 3 e 3 F 3 (r 3 ,y,e 3 ,p 3 ), (6.15) 

V = 1 — /3(1 — M3 )w(y) + 0(r 3 ), 
e 3 = -2e 2 3 F 3 (r 3 ,y,e 3 ,p 3 ), 

M3 = -e3P3F 3 (r 3 ,y,e 3 ,p 3 ), 

where we have divided the vector-field by r 3 and set F 3 (r 3 , y, e 3 , M 2 ) = 6 + aw(y) + 0{r 3 ). The analysis of 
the dynamics in this chart is very similar to the analysis in chart k\. Therefore we state the main results 
only. 
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There exists a repelling center manifold M r> 3 of r 3 = 0, y = w 1 , £3 = 0 given by 

M r ,3 ■ V = w _1 ( -577-^-c ) + 0 (r 3 + e 3 ) 

V/5(l — Ms) 7 

and a center manifold C rj 3 within r 3 = 0 which takes the form 

a ’ 3;S = ”"'(«rbi)) +0(es) - 

Within C rt3 there exists an invariant line 


lr ,3 : r 3 = 0 , y = y c , e 3 > 0 , y 3 = 0 . 


The center manifold C r)3 is overflowing and hence unique near l r 3 for II > 0. If Q < 0 then C r 3 is non-unique 
near Z 0j3 . The sub-manifold M ri3 (e) = M r>3 D {r 3 e 3 = e} is the extension of S r , f into chart n 3 near the line 
q. Also M rj3 (e) intersects 

A3“ : e 3 = i/ 2 , ( 6 . 16 ) 


which corresponds to Aj through the coordinate change K2 3 in (16.61) . ©(y^-close to the intersection of C r>3 
with (16.161) 


C r , 3 n A^ : e 3 = v 2 , r 3 = 0 ,y = y c + 0 (y 2 + y 3 ). 


The line l r<3 and the section A^" are shown in Fig. 16.41 

6.4. Persistence of canards. In this section, we combine results from the three charts «i _ 3 to show 
how singular canards in our PWS system can survive the regularization. For singular canards, there were 
two cases to consider: II > 0 (vrai singular canards) and < 0 (faux singular canards). We now consider 
the same two cases for canards in the regularized system. 
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For n > 0, C a , 1 and C r , 3 are unique as center manifolds. By transforming I 2 D of (16.121) into charts 
K\ and K 3 we conclude that K12U2) = la, 1 C C a . 1 for X2 <C 0 and K32{h) = lr, 3 C C r , 3 for X2 0. It is 
therefore possible to extend C a , 1 and C r , 3 into «2 as invariant manifolds C a , 2 and C r , 2 by using the invariant 
line I2 as a guide. Then from the analyses in section 16.21 and section 16.31 we conclude that the manifolds 
M a 2 (e) = K 2 i(M aj i(e)) and M r ^(e) = K 23 (M r> 3 (e)) intersect a distance O(r 2 )-close to the intersections 
of C a ,2 and C r> 2, respectively, where r 2 = yfe by (16.21) . By applying regular perturbation theory in chart K2 
the manifolds M at 2 (e) and M rt 2 (e) can be continued all the way up to *2 = 0 where they remain 0 (r 2 )-close 
to C a ,2 and C r , 2 , respectively. 

The manifolds C a ,2 and C r> 2 intersect along I2. We now investigate whether the intersection is transverse. 
If so, we can conclude that M a , 2 (e) and M r , 2 (e) also intersect transversally 0 (r 2 )-close to I2, since M a , 2 (e) 
and M r , 2 (e) are O(r 2 )-close to C a ,2 and CV-, 2 , respectively. We first eliminate time from (16.81) by division with 
X 2 and re-writting our equations in terms of w instead of y using (|4. 121) . Then we let (u{x 2 ), y- 2 ( 2 : 2 ), ^ 2 (^ 2 )) 
denote the variations of (w, ^2,^2) about I2 ■ 

X(/3X2'U — jx 2 

0 , 

0 , 

where 

A = G 2( x 2) = P- 1 (xl(v~ + PV + ) ~ Vc{x~ + Px + )) , (6-18) 

and 


du 
dx 2 
dfi 2 
dx 2 
dr 2 
dx 2 


= (/3 1 ,0,0). This gives 

+ f 2 G2{x2,u)), ( 6 - 17 ) 


4>i,c = ftiVc) > 0 . 

For simplicity we now drop the tildes from y 2 and r' 2 - We will apply the following lemma to (16.171) . 
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Lemma 6.1. (13i\ Proposition 4-%]) Ca ,2 and, C Ty 2 intersect transversally along I 2 if and only if there 
exists no non-zero solution of {6. 17 p | r2 — p which has algebraic growth for both X 2 ±00. 

Proof. Variations within the tangent spaces ToCa ,2 and TqC t ^ are characterized by algebraic growth in 
the past ( X 2 —> — 00) and in the future (a ;2 —> 00), respectively. Since C a ,2 and C r ,2 are unique, variations 
normal to ToCa ,2 and TqC t ^ will be characterized by exponential growth in the past and in the future, 
respectively. □ 

The following lemma describes the properties of the solutions of (IG.lTD r- _n necessary to invoke Lemma lG. 1 1 

Lemma 6.2. If Q > 0 then 0 has two linearly independent solutions 



The solutions 7 ^ grow exponentially as X 2 —► ±00 but the growth is only algebraic as X 2 —> Too, respectively. 
If II < 0, then neither of the solutions of {6.1 T^ -p grows exponentially as X 2 —> ± 00 . 

Proof. For II > 0 the result follows from the asymptotic behaviour of the error-function: 

9 r z 

erf(z) = —= / e _t dt. 

V 71 " Jo 

For It < 0 we use that f( e^^* 2 dt = 0(z~ 1 ) for z -A ± 00 . □ 

Therefore by Lemma Tb. II the manifolds C a ,2 and C r ,2 intersect transversally for H > 0 along ? 2 - Hence 
M a> 2 (e) and M r ^(e) intersect transversally (^(y^-close to r 2 for p, 2 = H 2 ,c where 


P 2 ,c = 0{r 2 ) (6.19) 

and r 2 = y/e sufficiently small. In fact, the following lemma allow us to calculate p, 2 ,c to lowest order: 

Lemma 6.3. Consider H > 0, then the only solution of {6.1 71 ) which has algebraic growth for X 2 —^ ±00 
is 


u{x 2 ) = -r 2 /3 2 (r] +/3ri + )x 2 , 
for 

M 2 = - (t^I + Pri + )+yc{x~ + /?X + )) /J~ 1 r 2 - (6.20) 

Proof. The complete solution of (16.171) is 

u(x 2 ) = ^ - (r 2 /? -2 (rj~ + N + ) + Kt 2 + r 2 fJ~ 1 yc{x~ + /3x + ))) \J^ eri + 

- r 2 fl~ 2 {r]~ + f3r/ + )x2. ( 6 . 21 ) 

Setting it(0) = 0 and (16.201) the result follows from Lemma RT21 □ 

Then following [20, Prop. 3.1] and using (16.201) we obtain 

M 2 ,c = ~ ( (^+^ 2 ^ C 7 ?" + P J 1 + ) + MX~ + Px + )^j P~ X r 2 + 0(r 2 2 ), (6.22) 


and hence we have the following main theorem. 
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Theorem 6.4. Assuming m , the regularized system m has a maximal canard at 

M = = °( e )> 

where /x 2iC is given by 16.22]) . which is 0(yfe)-close to the singular canard. 

Proof. The statements follow from the analysis above. □ 

In the case where 11 < 0 then C a ,i and C r ,3 are both non-unique and there exists a faux canard. The 
proof of this is very similar to the proof of faux canards in R 3 mm- 

7. Equilibria of the regularized system. We now revert to describing the dynamics in terms of the 
fast time r in (14.1511 for the remainder of the paper. In section 13.21 pseudo-equilibria of the PWS system 
were shown to exist for aS < 0 only. In Theorem l4.21 we showed that the equations of the reduced problem 
(14.71) agree with the sliding equations (12.141) and (12.151) . Hence, as a consequence of the implicit function 
theorem, pseudo-equilibria of the sliding vector field X s ;(x, p) for \p\ > po > 0 perturb to locally unique 
equilibria on the slow manifold of the regularized problem for 0 < e < eo(^o) sufficiently small. Indeed, from 
our fast system (14.161) . we find an equilibrium at 

x* = pSfl-'n + Oin 2 +e), ( 7 . 1 ) 

y* = w~ 1 (—a~ 1 5) + 0(p + e). 


Since we require w(y) > 0 by (14.121) . (x*,y*) is an equilibrium provided aS < 0. This is in accordance 
with equation (13.21) for the PWS system. The stability of this equilibrium is described by Proposition 13.31 
However, for p = 0, e = 0, the equilibrium lies on the non-hyperbolic line q , defined in (14.111) . of the critical 
manifold. Hence Fenichel theory can not give a description of the equilibrium for e ^ 0. To accurately follow 
the equilibria near p = 0 we need to consider the extended equations in chart «; 2 , equations m ■ 

7.1. Chart k 2 . From (16.81) . we find the following equilibrium: 


(x%,y*) = (P 6 Q. 1 H2,w t (-a 1 S))+0(r 2 ). 
The equilibrium (a" 2 ,y*) intersects the section Aj, defined in (16.111) . when 

y* = w _1 (—a _1 <5) + 0 (r 2 ), 

H 2 = + 0 (r 2 ). 

Consider r 2 = 0. The linearization of (16.81) about (x 2 ,y*) is then given by 


z = Az, A 


1 ) aw'(—a 3 (5) 

a/3fl -1 u/(— a~ 1 S)y.2J ’ 


where z is the variation of (x 2 ,y) about (a; 2 ,y*) and 

mr) 


w'(—a M) = — 
using (14.131) and the fact that 


1 




(7.2) 

(7.3) 


(7.4) 




1 + a 

1 — a~ 1 5 


The determinant of the coefficient matrix A is independent of /r 2 : 

det A = — Clw'(—a~ l 8 ) 


and the trace of A is given by 


trA = a/3fl 1 w'(— a 1 5)/Lt 2 


(7.5) 

(7.6) 


which vanishes for /i 2 =0. Since w' < 0 from (14.131) . the sign of det A is determined by the sign of fi. 
Therefore we conclude the following: 

Lemma 7.1. Consider r 2 = 0 and suppose aS < 0. 
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• For Q < 0 we have: 

— The equilibrium (x 2 ,y*), defined in is a neutral saddle for p 2 = /j. 2 ,n = 0 and a generic 

saddle for y, 2 ^ y .2 ,n ■ 

— The stable (unstable) eigenspace associated with (x 2 l y*) and the linearization ra is asymp- 

totically vertical for —>■ Too. 

— The unstable (stable) eigenspace associated with (x 2 l y*) and the linearization (Z2P is asymp- 
totically horizontal for /3<5S4 _1 /r2 —> Too. 

• For Q > 0 we have: 

— A Hopf bifurcation at y 2 = ^2 ,h = 0. 

— The equilibrium (x 2 ,y*) is attracting (repelling) for (dSTL -1 y .2 ;s 0. 

— For \y 2 \ e (V 2 ,H,ti 2 ,F), where 


= 20 3 / 2 
/i2 ' F “ \J—w'(—a- 1 8)\a!3\ ’ 

t/ie equilibrium (x 2 ,y*) is a focus. 

— For |//21 > H 2 ,f, {x 2 , y*) is a node. 

— The strong eigenspace is asymptotically vertical for /3(5fi _1 /Lt2 —► ± 00 . 

— T7ie wea/c eigenspace is asymptotically horizontal for (55Tl~ l y 2 ±00 

The quantities y. 2 ,N> M2 ,h and y 2 ,f perturb by an amount of 0(r 2 ) for r 2 7^ 0 by the implicit function 
theorem. 

Proof. From (17.411 we find the following eigenvalues and eigenvectors: 

A± = ^ \/tr 2 A — 4det A, (7.8) 

v± = (f) ' (™) 

For p .2 = 0 we have a neutral saddle for B < 0 and a center for fl > 0, since tr A = 0 from (17.(ill and the sign 
of det A is determined by the sign of 0. Also, since T^-tr A 0, we have a Hopf bifurcation at M2 = P 2 ,h = 0 
when fl > 0. 

For fi > 0 we note that when 


tr A = ±2Vdet A, 

that is, when y 2 = TM2 ,f from (ED , then the equilibrium (x 2 ,y*) is a node. It is attracting (repelling) for 
those fj .2 for which tr A ^ 0. So from (17.61) . we conclude that {x* 2 ,y*) is attracting (repelling) for (3SQ 1 y 2 0. 
For either sign of H we have, for /3<5f2 -1 /i2 —> ±00, that 

A + —y Too, 

A_ -r 0. 

So, for S2 < 0, using urn this means that the stable (unstable) eigenspace associated with (x 2 ,y*) and the 
linearization m is asymptotically vertical for (35Q 1 y .2 —> Too- On the other hand the unstable (stable) 
eigenspace is asymptotically horizontal for /35H _1 /i2 —>• Too. 

Similarly for H > 0 we find that the strong (weak) eigenspaceQ associated with (x 2 ,y*) and the lineariza¬ 
tion ED is asymptotically vertical (horizontal ) for /3<5f2 1 y 2 —i * ±00. □ 

The results of Lemma T7. II are shown in Fig. 17.11 which can be compared with Fig. 13.21 

7.2. Charts ^1,3. The results in Lemma 17.11 are in accordance with Proposition 13.31 in the limits 
y .2 —> ±00. But they occur within chart k 2 and everything collapses to y = 0 for r2 = 0, see m - To 
connect the results in chart with the case fi = 0(1), we can consider charts We obtain the following: 


6 The strong (weak) eigenspace is the eigenspace associated with the eigenvalue representing the stronger (weaker) contraction 
or expansion. 
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/35ft 1 H2 -C 0 


ft < 0 



ft > 0 





v = y* 


^2 


/35ft 1 fi2 — 0 


ft < 0 



X2 


ft > 0 



y = v 


X2 


/35ft 1 /x 2 » 0 


ft < 0 


ft > 0 




Figure 7.1. Results of Lemma \7. 1\ The dashed line corresponds to y* = w 1 (—o: 1 6 /, from f 7. 2| ) . T/iis figure can be 
compared with Fig. 13.21 


Lemma 7.2. TTie equilibria described in Lemma \7.1\ belong to the same smooth, locally unique family of 
equilibria as those in \7.1\ . The equilibria do not undergo any further bifurcations in passing from A 2 to 
A ± . 

Proof. In chart n \, we find the following family of equilibria: 

y = ur 1 (-a~ 1 S) + <D(n + a), 

Mi = 1 + 0 (r\ + ei), 

within M a \. Using the conservation of e and /i it is straightforward to trace this family of equilibria from 
Af to A - and connect them to the equilibria described in chart with the 0(1) equilibria described in 
(O). The analysis in chart K 3 is identical. □ 

7.3. The Hopf bifurcation for fl > 0. We now describe in further detail the Hopf bifurcation for 
ft > 0 in Lemma P and the resulting birth of limit cycles. As mentioned at the start of this section, we 
return to the time r, since time t defined in (14.151) and used in section [6] leads to difficulties when the periodic 
orbits leave the region of regularization y S (—e, e), or y £ (—1,1) (whereas canards do not suffer this fate). 
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To proceed, we write the extended problem (14.1411 in chart k 2 to obtain 

x 2 = (6 + r 2 ( + x 2 + 0(r 2 (p 2 + r 2 )))(l + 0{y)) + (a + r 2 C,~x 2 + 0(r 2 (p 2 + r 2 )))(l - 0(y)), (7.10) 

V = {x 2 + r 2 r] + x 2 + 0(r 2 (y 2 + r 2 ))))(l + 0{y)) 

+ (-/?( x 2 - fj, 2 ) + r 2 r]~(x 2 - M 2) 2 + 0(r 2 (n 2 +r 2 )))(l - 0(y). 

where we have multiplied time r by r 2 . The equilibrium in (17.21) for a.5 < 0 then becomes 

(z 2 ,y*)= (j )- 1 + 0(r 2 ). ( 7 . 11 ) 

This equilibrium undergoes a Hopf bifurcation at y 2 = p 2 ,H = 0(r 2 ) when Q > 0 (compare with Lemma l7.ll) . 
Let 


where 


and 


0i,H = 0 M (yo), 


* = 1,2,3, 


(7.12) 



/l + a-M\ 

Vl-a-My ’ 


(7.13) 


d l 0 


The subscript 0 in (17.131) is used to emphasize that y* has been obtained from (17.111) with y, 2 = r 2 = 0. By 
assumption <pj t H > 0 and so we obtain 

Proposition 7.3. System 1 7. 1 U\ ) has an equilibrium at (x 2 ,y*), as defined in 1 7.11\ ), which undergoes 
a Hopf bifurcation at 


M2 = M2 ,h(x 2 ) 


where 


M2,ff(T 2 ) 


/ 2 ( a (c+ + x + )-«y(C~+x~))» 

v {a - 5) 2 (fiM 


(x + Px Wo )P r 2 + 0{r 2 ). 


The first Lyapunov coefficient is given by 


(7.14) 


a = a 2 r 2 + 0(r 2 ), 


where 


a 2 = 


+ 


(a - 6)0\ H 


16S2 2 
1 I 03, H 


— ^(^ + 1) 2 (^C ~ a( + ) - {a - 5) 2 (p + Py + ) ^ ^ + iq^’ H (<*X ^ a X + ) 


b 2 ,H 


8 l 0i, h( a - 6) 0\ H (a - 6) 


(/? + ^ 2 ’ g ) ( S(C + X_) “ a(C+ + ^ + )) • ^ 7 - 15 ) 


7/ a2 / 0 then for e sufficiently small there exists a family of unique periodic solutions bifurcating from 
{x 2 ,V*) for 

fdSVT 1 p 2 > f35Q~~ M2,£f when a 2 < 0, 

/3<5fl _1 /i 2 < P6fl^ 1 p 2 ,h when a 2 > 0, 

with amplitude O (^J\p 2 — M 2 ,i/|t^ 1 ^ • The periodic orbits are attracting for a 2 < 0 and repelling for a 2 > 0, 
for r 2 sufficiently small. 
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Proof. The calculation of 02 is based on classical Hopf bifurcation theory [5|. □ 

Note that the Hopf bifurcation is degenerate within r 2 = 0 since a = 0 there. The reason for this is that 
the system (17.101) with p 2 = 0 and = 0 is Hamiltonian, as we shall now demonstrate. In this case, (17.101) 
becomes: 


x 2 = 5(1 + <j>{y)) + a(l - 0(y)), 

y = (1 + 0(y) - /3(i - 0(y))) * 2 , 

The Hamiltonian function H = H(x 2 ,y) is given by 

5(1 + 0(s)) + a(l - 0(s)) 


1 f y 

H(x 2 , y) = -x\ + / 
z JvZ 


'vs - ( 1 + 0 ( s )) 

The symplectic structure matrix J(x 2 ,y) is non-canonical: 


ds. 


J(x 2 , y) = (/3(l - 0(j/)) - (l + 0(y))) 


0 1 

-1 0 


(7.16) 


(7.17) 


(7.18) 


which is regular and non-zero near y = y g since 0 (j/q) = , aS < 0 and Assumption [3] Hence the 

system with /i 2 = 0 and r 2 = 0 has a whole family of periodic orbits in the vicinity of (17. 1 ID ...___n. The 

Hamiltonian system is not well-defined for y = y c (|6.10l) . /3 > 0, since J{x 2 ,y c ) = 0. 

Remark 7.4. The periodic orbits within the (x 2 ,y)-plane rotate about \1.11 1 1 in the counter clockwise 
(clockwise) direction if a > 0 (a < 0 ). 

Combining the results in Lemma l7.ll Lemma l7.21 and Proposition [773] we obtain one of our main results, 
Theorem 17.51 as follows: 

Theorem 7.5. Assuming \5.2\) . the regularized system j4-3\ ) h as a smooth and locally unique family of 
equilibria 


{x, y) = (x*,y*)(p, \fe) = V + + e), O(e)) 

where (x*(y,0),y*(y 1 0)) = (x*(p.,0),0) agrees with the family of pseudo-equilibria for the PWS system (see 
Provosition \3.3\ ) and, in particular, d^x*(p,0) 0. For 

fl < 0: The family of equilibria consists of saddles and does not undergo any bifurcation. 

12 > 0: The family of equilibria undergoes a Hopf bifurcation at p = \fep 2 ,h(\/c) = 0(e), where p 2 ,h is 
given by \7.ll$ . The first Lyapunov coefficient is given by a = a 2 \fe + 0(e) where 02 , given by 
1 7. 15\ ), depends upon the regularization function <j>. 

Note that 02 , as given in (17.151) . depends upon the regularization function <p, through <f>i } H, 4>2 ,h and f> 3 t H 
as defined in (17.121) . and hence that the criticality of the Hopf bifurcation depends on <f>. This observation 
leads to another one of our main results: 

Theorem 7.6. Suppose that Q > 0 and aS < 0 so that there exists a Hopf bifurcation. Then, provided 

5(C+X”)-a(C + +X + )^0, (7.19) 

all cases <22 < 0, 02 = 0 and 02 > 0 can be attained by varying (f>. 

Proof. We are free to choose <j>\> 0, (j> 2 ,h and in (17.151) . Note, from (17.151) . that the equation 
02 = 0 is linear in f> 3 ,H and that the coefficient of 

5(C~ +X~) ~ «(C + +X + ) 

80i ,h{ol - 5) 

is non-zero by assumption (17.191) . Hence 02 = 0 can be solved for 03, h- The statement of the proposition 
therefore follows. □ 

Remark 7.7. It seems natural to insist that 0 should be an odd function. If 0 were not odd, then one of 
the vector-fields would be favoured over the other by the regularization. The functions 0 that we used to 
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prove Theorem | 7. 6 \ can be odd, at least if y^ 7 J 0. J/yJ = 0 and <j> is odd, then 4>2,h = 0. Hence, the equation 
a 2 = 0 should have a solution with <p 3 ,H < 0 for <f> to be odd and for Theorem m to apply. 

Remark 7.8. Another natural condition appears to be that < f >' should be strictly increasing within (—1,0) 
and strictly decreasing within (0,1). The functions used in the proof of Theorem | 7. (i\ may also be chosen to 
satisfy these conditions, at least when y J 7 ^ 0. The functions then just have 4 > 2 ,h S® 0 for y* ^ 0. For y* = 0 
we have 4 > 2 ,h = 0 and the equation 02 = 0 should have a solution with 4 > 3 ,h < 0 to ensure 4 >"{y) ^ 0 for 
y gs 0 and that Theorem | 7. 6 '| apply. 

8. Limit cycles of the regularized system. From section ITTTTT1 limit cycles of the PWS system can 
exist for the case II 2 , which occurs when 

S = — 1, a > 0, /3 < 0. (8.1) 

See also Table O Limit cycles can also occur in the regularized version of case VI 3 . This case occurs for 

5 = 1, /? > —a > 0. (8.2) 

Note that /3 > 0 in (18.211 and hence from Theorem [631 the regularization of VI 3 also possesses a canard which 
is 0(v / e)-close to the singular canard of the PWS system. 

The regularization of II 2 and VI 3 , denoted by JJJ and VJJ respectively, exhibit significant differences. 
In case Ilf, the limit cycles eventually (for y 2 large enough) cross the region of regularization y € (— 1 , 1 ) 
from y > 1 to y < — 1 and back again. There is no sliding and hence no singular canards in the corresponding 
PWS case I/ 2 . Thus there are no canards in the regularized case JJJ. However, for VJJ, the resulting limit 
cycles interact with the slow manifolds and the maximal canard to produce a scenario almost identical to 
the canard explosion phenomenon in classical slow-fast theory [ 21 ] . 

In this section we present a comprehensive study of the regularized limit cycles that are due to the 
Hopf bifurcation in chart k 2 (see Proposition 17.31) . In section fTTl these limit cycles are followed, beyond 
the validity of the classical Hopf bifurcation theory, into large limit cycles in chart K 2 - In terms of the 
original (x, y)-variables, from (16.21) . these periodic orbits are, however, still small, only extending 0 (^/e) in 
the cc-direction and O(e) in the y-direction. To follow these orbits to (H(l)-size, and obtain a connection to 
the PWS system, we must use charts /ci^. 

In doing so, we use different techniques for cases JJJ and VJJ. We split the analysis into separate parts. 
In section l8.2.1l we study limit cycles of CJ(l)-size for case JJJ while section 18.2.21 contains the corresponding 
analysis for case VJJ. The connection of these 0(l)-limit cycles with the limit cycles in chart n 2 for cases 
//J and VJJ is shown in sections [8.3.11 and 18.3.21 respectively. 

8.1. Chart k 2 . Proposition 17.31 only guarantees the existence of small periodic orbits for small 

within chart n 2 - To follow these periodic orbits within chart n 2 for larger values of ^ 2/^2 we follow the 
Melnikov-based approach of Krupa and Szmolyan m • We will consider both J/J and VJJ in this section. 

First we define yj, y f and y 2 as follows. Consider the forward solution a = (x 2 ,y){t) with initial 
condition (0, y§) where the Hamiltonian H in (17.171) takes the value JJ(0,yo) = h > 0 and 

Vo<Vo- (8-3) 

The point (0, y 2 ) is then the second return of er to x 2 = 0 where 

y -2 < Vo- 

Similarly, we let (0, y^) denote the first return of er to x 2 = 0 where 

Vi > Vo- 

The relevant quantities are shown in Fig. 18.11 for the case VJJ. Notice that by (18.21) and ST > 0 it follows 
from (16.91) and (17.21) that the singular canard for y 2 = r 2 = 0 is above the bifurcating equilibrium. This is 
illustrated in Fig. 18.11 bv letting the continuation of S a ,e and S r>e lie above yj. The case JJJ is similar but 
there are no slow manifolds in this case. 
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Figure 8.1. Illustration of y(j, y’f and ylf for the case V/J. Here S ai e and S r ,e are in fact the continuation M a , 2 (e) and 
M r ■. 2 (e) of the Fenichel slow manifolds into chart K 2 - The case //| is similar but there are no slow manifolds in this case. 


Following Proposition 17.31 there exists an Jig > 0 independent of r 2 so that for 0 < h < ho there exists a 
locally unique family of limit cycles parametrized by /i2 whose stability is determined by the sign of 02 given 
in (17.1 51) . We therefore take h > ho > 0 and consider the following distance function 

D(r 2 ,fi 2, h) = H{0 , y£) - H( 0, y h 0 ). (8.4) 

From the analysis proceeding Proposition 17.31 the system with y ,2 = ^2 = 0 is Hamiltonian and hence 
D{ 0,0, h) = 0 for all h > 0. Also since 

* H = ^ + </>($)) + a 0 - </>((/)) , n 
" /3(1 — <j>{y)) — (1 + * ’ 

for y < Vo (in accordance with (18.311 ), roots of the equation D = 0 correspond to periodic orbits. 

Let T h denote the period of the orbit a of the Hamiltonian system with /i2 = 7~2 = 0 satisfying H(a) = h. 
We have the following lemma, similar to [211 Proposition 4.1]: 

Lemma 8.1. For h > ho we have 

D(r 2 , /i2, h) = D T 2 (h)r 2 + 2 (h)/i 2 + 0((r 2 + h 2 ) 2 ), (8-5) 


where 


D r2 (h) = 2/ (C + (i + &{y)) + C (1 

JO 

S(1 + <p(y)) + a(l ~ 


xl 


+ 


= -4 


PO - <£(£)) - (1 + 4>{y)) 

W{y) 


((v + xl + x + y)0 + <f>(y)) + (v %t + x y)0 - 4>{y)))dt, 

- X 2 dy, 


( 8 . 6 ) 


’Vi 


Dfj, 2 (h) 


m-m)-o+m )) 2 


(8.7) 
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where (x 2 ,y)(t) satisfies \7.16\ ) with (x2,y)(0) = (0 ,$), H(x 2 ,y) = h and (x 2 ,y)(T h / 2 ) = (0,S/i)• 
Proof. Similar calculations to those in ED Proposition 4.1] lead to 


D r 2 {h)= I VH(x 2 (t),y(t)) ■ G r 2 (x 2 (t),y{t))dt, 
Jo 

H>uA h ) = [ VH{x 2 (t),y{t)) ■ G^{x 2 {t),y(tj)dt, 

Jo 


where 


Gr 2 (X 2 ,y) = 
G/xs i x 2,y) = 


(C + (i + 4>{y)) + C (i - <t>(y)))x 2 
(y + xl + x + y)(l +4>{y)) + (tT®! + x~ti)0 ~ Hi)))J ’ 

° ) 

(1 - <f>{y))fi) 


The Hamiltonian system possesses a time-reversible symmetry (x 2 ,y,t) K > (—X 2 ,y, —t) so: 

f T h /2 

io 


D r 2 (h) = 2 / VH(x 2 {t),y(t)) ■ G r 2 (x 2 {t),y(t))dt, 

Jo 

Dfi 2 (h) = 2 / S7H(x 2 {t),y(t)) ■ G fl 2 (x 2 (t),y(t))dt. 

Jo 


For U M2 (h) we then use integration by parts. □ 

Remark 8.2. If \yo\,\yi \ > 1 then can be simplified further: 


D„Jh) = -4 J 


fix 2 


-i 




( 8 . 8 ) 


This is only relevant for case Ilf,. In the case Vlf the maximal canard prevents the local limit cycles from 
entering y > 1 (see Fig. EIp. 

Since 2 (/i) ^ 0 we can apply the implicit function theorem to conclude the following: 

Proposition 8.3. Fix v > 0 small. The family of limit cycles from the Hopf bifurcation, described in 
Proposition | 7. S\ can be continued into periodic orbits corresponding to roots of D(j 2 , P 2 , h) fork < 2v~ 2 and 
r 2 < r 2n(p) sufficiently small. The orbits are 0 (r 2 )-close to the periodic orbits of the Hamiltonian system 
H = H(x 2 ,y) defined in [7J7\) . 

Proof. From (18.51) and the implicit function theorem we obtain 

»* = -w(M r2+0(rl) ' (8 - 9) 

□ 

Remark 8.4. The periodic orbit with h = 2z/~ 2 intersects y = j/q in x 2 = = f2i /_1 for P2 = r 2 = 0. By 
Proposition ^. fH we are therefore able to continue periodic orbits beyond the sections AJ in chart H 2 - These 
orbits belong to a C k -smooth and locally unique family because they are obtained using an argument based on 
the implicit function theorem for y .2 = ^2 = 0. 

8.2. 0(1) limit cycles. We now wish to consider limit cycles of the regularized system with amplitudes 
that are 0(1) with respect to e. As mentioned above, the analysis is divided into two cases: Ilf and V7J. 

8.2.1. 0(l)-limit cycles for case Ilf. We start by obtaining 0(l)-periodic orbits in the original 
(a:, y)-variables by following fixed points of a Poincare map: 


P e : {y = e, x > 0} —)> {y = e, x > 0}, 


( 8 . 10 ) 


where defined under the flow of the regularized system. Since these orbits are 0(1) with respect to e and 
only involve crossing (see Fig. 13.31 case II 2 ), the mapping P e is smoothly 0(e)-close to Pq, as defined in (13.61) 
for the PWS system. We therefore obtain the following proposition. 
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Figure 8.2. Important quantities for case V/J, relevant to Proposition [8f6\ The region of regularization y E (—e,e) is 
exaggerated for illustrative purposes. 

Proposition 8.5. Fix n$ small and consider m - Suppose that |/z| > po > 0 and A jj 2 ^ 0, where 
A jj 2 is defined in \3.8\) . Then for e < eo(po) sufficiently small, the regularized system has a family of periodic 
orbits corresponding to fixed points of P e of the following form 

x e (p) = x 0 {p) + 0 (e), 

with Xo(p) given by S3.9\) . and pAj^ < — /xo|A// 2 | _1 < 0. The periodic orbits are attracting for A j / 2 < 0 and 
repelling for A jj 2 > 0. Moreover, they are continuously 0(e)-close to periodic orbits of the PWS system. 

Proof. Proposition 18.51 gives non-degenerate fixed points Xq( p) of Po for pAj^ < — p 0 \ A// 2 | 1 < 0 . 
Since P e = Po + 0(e) we can apply the implicit function theorem to obtain fixed points x e (p) = xq( p) + 0(e) 
of P e . The stability of the fixed points x e (p) is also determined by the stability of xo(p ) as a fixed point of 

Po- □ 

8.2.2. 0(l)-limit cycles for case V7J. This case has a canard at p = r 2 P 2 ,c and a Hopf bifurcation at 
p = r 2 P 2 ,H■ The co-existence of a Hopf bifurcation and a canard leads to the canard explosion phenomenon 
in which the amplitude of limit cycles undergo 0 ( 1 ) variations within an exponentially small parameter 
regime. In order to prove this statement, we follow the proof of a related assertion in [ 2T| Proposition 5.1] 
for classical planar slow-fast systems. 

Let p = r 2 P 2 = \ftP 2 and consider the original (x, y)-variables, in which the region of regularization is 
y € (—e, e), and let 7 be the forward orbit with initial condition (x,y)( 0 ) = ( 0 , yo) where yo < 0 small but 
independent of e. The situation is illustrated in Fig. 18.21 Since the fold of the associated PWS system is 
invisible from below (see (18.11) and Proposition 12.41) we know that the first return of 7 with x = 0 is a point 
(0,yi) with 2/1 = 0(e)- Let 7 be the backward orbit with initial condition (x, 2 /)( 0 ) = (0, j/o)■ Denote by 
(0,2/ 1 ) the first return of 7 with x = 0. Here y 1 = 0(e). Let (xi(?/o, e), —e) and (^ 1 ( 2 / 0 , e), — e) denote the 
first intersections of 7 and 7 , respectively, with y = —e. The functions xi( 2 /o, e) and x^o, e) are smooth in 
2 /o and e. In particular, xi( 2 / 0 , 0 ) and x^j/OjO) can be obtained from the associated PWS system. 

We consider the distance function: 


V(y 0 ,p 2 ) = 2/1 - V v 

Roots of V correspond to periodic orbits. As in m we solve T>(yo , P 2 ) = 0 by noting, from Fenichel theory, 
that 


r>{yo, p 2 ,c) = 0(e -c/e ) 


c = c(//o) > 0. 
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Since SV,e and S a>e are transverse for y 2 = y 2 ,c this then effectively implies the existence of a /12 = M 2 ,p(j/o) 
solving U(yo,[i 2 ) = 0 and satisfying y 2 , p = y 2 ,c + 0(e~ c ^ e ). Since X\ and are increasing functions of 
yo for ?/o small, it also follows that y 2 , P (yo ) approaches /i 2 ,c monotonically as yo increases, at least for 7/0 
sufficiently small. The stability of the periodic orbits is determined by the sign of a way-in/way-out function 
R = R(yo ), see |21[ Proposition 5.4], that measures the contraction and expansion along S atr . In our case 
the contraction and expansion is determined by the following function 

(X£{x, 0, 0) — Xf (x, 0, 0))4>'(y). 

Inserting y = y c = 1 j)^ 1 from (16.101) . which corresponds to S a , r for /i 2 = £2 = 0, we obtain the function 

R = R{y 0 ): 

r2.i(yo,0) 

R{yo) = / (-^ 2 ”( x > 0,0) — Xf (x, 0,0))d£. 

dxi(i/o,0) 

Since <j> i jC = <f>\y c ) > 0, the sign of R coincides with the sign of 

rn 1 ( 2 / 0 . 0 ) 

R(y 0 ) = / (X+(x, 0, 0) - Xf(x, 0,0 ))dx. 

Jxi(y 0 ,0) 

We obtain the following proposition: 

Proposition 8.6. Consider a point p = (0, yo) with yo € [—c]" 1 , — cf 1 ] with C 2 > c\ sufficiently large but 
fixed. Then for e < eo(ci, C 2 ) sufficiently small there exists a unique periodic orbit through p for y = y p (r 2 ,yo) 
where 


H P (r 2 ,yo) = he + 0 (e c/e ). 


The function y p (r 2 ,-) monotonic so that y p (r 2 ,yo) approaches y c {r 2 ) as yo increases. 

The periodic orbits are attracting if 

Avh = 3af3 + ^ V+>> + ^ + + ’ ( 8 - n ) 

is negative. They are repelling if A yi 3 is positive. 

Proof. To verify the statement about stability we need to compute R(yo). To do this we invert X\ = 
£ 1 ( 7 / 0 ,0) for yo and parametrize x x in terms of £1 rather than 7 / 0 . Then x^xi) is obtained from the map afi 
in Lemma 13.51 with y = 0: 

£1 = crfi{xi) = —£1 + A~ x\ + 0{x\), 


using the backward flow of X , where A is given by (E3]). We therefore consider the following integral 

rx 1(2:1) 


R{x 1 ) = 


(X+(£,0,0) - A 2 "(£,0,0))d£ 


rx i(zi) 

/ ((1 + (3)x + (rj + — y~)x 2 + 0(x 3 ))dx 

J 


= ( A (1 + /3) + — {rj + — 77 ) |£i | 3 + (D(xj). 


Hence the sign of R{x 1 ) is determined by 


A (i + fi) + ^(v + ~ v ) = (a(y + Pn + ) + P{/3 + i)C ) 


where we have used (13.51) . The right hand side is (18.111) . Since R < 0 implies stability while R > 0 implies 
instability the result follows for X\ (and hence yo) sufficiently small. □ 

Remark 8.7. In classical planar slow-fast systems m, a canard is generically associated with a Hopf 
bifurcation and a canard explosion. This is not necessarily the case here. For example the PWS case VV\ 
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has a singular canard but no local limit cycles in either the PWS system or its regularization. Conversely, 
the regularized system can undergo a Hopf bifurcation without the presence of a canard. This is demonstrated 
by case J/|. 

We conclude this subsection with Table IO which summarizes properties of the seven regularized two¬ 
folds (compare with Table 13.111 . 


Table 8.1: 1: type of regularized two-fold singularity. 2: value of 
aS. 3: equilibrium (x = no, / = yes). 4: sign of fl. 5: Hopf (x = 
no, V = yes, n.a. = not applicable). 6 : value of j3. 7: canard (x = 
no, / = yes). 


1 

2 

3 

4 

5 

6 

7 

Type 

aS 

Equilibrium 

fi 

Hopf 

P 

Canard 

VVf 

+ 

X 

+ 

n.a. 

+ 

/ 

Wf 

- 

✓ 

- 

X 

- 

X 

VI{ 

+ 

X 

± 

n.a. 

- 

X 

VII 

- 

/ 

- 

X 

+ 

/ 

VII 

- 

✓ 

+ 

/ 

+ 

/ 

III 

+ 

X 

- 

n.a. 

+ 

/ 

III 

- 

/ 

+ 

/ 

- 

X 


8.3. Connecting limit cycles. Having obtained 0(1) limit cycles in the two cases 11% and VI%, we 
now analyze the connection between these limit cycles and those described by Proposition 18.31 that are due 
to the Hopf bifurcation. The results are summarized in one of our main results: 

Theorem 8 . 8 . For e sufficiently small: 

11%: There exists a C k -smooth family of locally unique periodic orbits of the regularized system \f.3\) that 
is due to the Hopf bifurcation in Theorem | 7. 5\ If 02 < 0 (a-i > 0,) where 02 is the first Lyapunov 
coefficient as defined in {7.15) , then the periodic orbits are attracting (repelling) near the Hopf 
bifurcation. If 

A// 2 = (oifq- + (3rj + ) + + x- + a(C + + X + ))) , (EH) 

is negative (positive) then the periodic orbits for A jfp < —c^fe, c > 0 sufficiently large, are attract¬ 
ing (repelling). The periodic orbits for fi = 0(1) (with respect to e) are continuously O(e)-close to 
the PWS periodic orbits in Provosition \3.6\ 

V /J: There exists a C k -smooth family of small periodic orbits of the regularized system )4-3) that is due 
to the Hopf bifurcation in Theorem \ 7.5\ There also exists a C k -smooth family of periodic orbits 
that are 0 ( 1 ) (with respect to e) in amplitude and which undergo a canard explosion, where the 
amplitude changes by 0(1) within an exponentially small parameter regime around the canard value 
H = \fep 2 ,cl\[e) (see Theorem W4\ >- If a 2 < 0 (a 2 > 0) where 02 is the first Lyapunov coefficient as 
defined in 1)7.15]) , then the periodic orbits are attracting (repelling) near the Hopf bifurcation. If 

Avh = 3o/3 + + M 3 + + )) d 811 D 

is negative (positive) then the O(\)~periodic orbits are attracting (repelling). 

The proof of Theorem 18.81 is divided into two cases: 11% and VI%. 

8.3.1. Connecting limit cycles for case II% using charts ^ 1 , 3 . To connect the periodic orbits 
described in chart K 2 by Proposition 18.31 with the 0(1) periodic orbits in Proposition 18.51 we first return to 
the original (a;, y) variables in which the region of regularization is y £ (—e, e) and consider the mappings of- 
taking (xo,e) to (cr+(cco),e) and (xi,— e) to (cr~(x 1 ),—e), respectively. See Fig. 18.31 We have: 
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Figure 8.3. The mappings erf associated with the regularization of the invisible two-fold singularity II 2 , where S = — 1 , 
a > 0 and (3 < 0. The points (x,y) = (rr^ibe) are points where X±(x, ±e, p) = 0. Note that the mappings and £ 3 " in 
( 8 . 16 \) and \8.11 p describe the fast transitions indicated in the figure by the triple-headed arrows, using the charts and K3. 
These transitions would appear nearly vertical in the stretched coordinate system ( x,y ). 


Lemma 8.9. Fix c > 0 large. The maps erf- are defined for 

x G (xy (yt, e), c~ l ] where x~j (/x, e) = — x + e + 0(e(p + e)), 
and 

x € [— c , xj (/z, e)) where xj (/x, e) = p — /3“x _ e + C>(e(p + e )X 
respectively and for those x the maps af satisfy 

=°o+ C’(e), 


( 8 . 12 ) 

(8.13) 


(8.14) 


where CTq are described in Lemma\3J^ 

Proof. The mappings cr± map {y = 0} to itself by the forward flow of X ± for x > 0 and x < p, 
respectively. The mappings cry, on the other hand, map {y = ±e} to itself by the forward flow of X , 
respectively. Here cr+ is defined for x > Xj(p, e) where X^ixpp, e), e, /x) = 0 while cr+ is defined for 
x < xj(p,e) where Xf(xf (p,e),—e, p) = 0. Using (12.12[l and (12.131) and the implicit function theorem gives 
(18.121) and (18.131) . respectively, for c sufficiently large. Equation (18.14f) therefore follows by standard regular 
perturbation theory. □ 

Remark 8.10. The mappings erf- do not depend upon the regularization. They are due to ra deter¬ 
mined by X ± . 

We now write these mappings af in terms of the charts k\ 3. The resulting mappings will be denoted 

by 


0-31 = k 3 o a e o Kl , 


and 


CT-^2 — Av]_ O CTg o 


respectively, using the subscripts to highlight that these mappings are from (k^) to K 3 («i), respectively. 






















36 


Lemma 8.11. Consider m- In terms of the charts k ±^ 3 the mappings af take the following forms: 

a 3 i '■ K 1 n {y = _1 } K 3 n {y = -1}, 

r 3 = n (1 + 2pi + A~n + 0(n(ei + pi +ri)), 

e 3 = ei(l — 2(2/ii + A n) + 0(ri(ei + Mi + r\) + p\)), 

M 3 = It l(1 — (2/^i + A r 3 ) + 0(r\(ei + Mi + Pi) + p\)), 

and 

vf 3 ■ « 3 n {y = 1} -A ki (~l {y = 1}, 

r\ = r 3 (l - A + r 3 + 0(r 3 (e 3 + n 3 + r 3 )), 
ei = e 3 (l + 2 H+r 3 + 0 {r 3 {e 3 + p 3 + r 3 ))), 
pi = p 3 {l + A + r 3 + 0(r 3 (e 3 + p 3 + r 3 ))), 
respectively, for rq , e\ and pi sufficiently small. 

Proof. Consider the case 0-7 (the case af is identical). We then use (18.141) and Lemma 13.51 and set 
x\ = — r 3 and x 3 = r 3 as described by the charts ki, 3 , respectively. The expressions for e 3 and p 3 then follow 
from the conservation of e and p, respectively. The condition (18.121) and (18.131) are satisfied for r 3 , e 3 , p 3 
and r 3 ei, pi, respectively, sufficiently small.D 

We then consider the Poincare mapping P e from (18.101) used in Proposition 18.51 and write this mapping 
in chart n 3 . The resulting mapping is given by 

P 3 = n 3 o P t o ■. k 3 D {y = 1} ->• k 3 n {y = 1}, (8.15) 

{r 3 ,e 3 ,p 3 ) (rt, 4 :ht), 

We will compose P 3 into four different mappings erjj, a 3 i and £ 3 " so that: 

P 3 = £+ o a 31 o o cr+. 

The mappings a 33 and are described in Lemma 18.111 while £7 and £ 3 " are given in terms of the forward 
flow associated with the differential equations in charts «i and n 3 (see (16.131) and (16.151) ') and map {y = ±e} 
({ y = ±1}) to {y = =Fe} ({y = =F 1}), respectively. Hence the mappings: 

Cf :{y = l}^{y = -l}, (8.16) 

(ri,ei,pi) (rf,ef,pf), 

and 

4 + : {y = -1} {y = 1}, (8-17) 

(r 3 ,e 3 ,p 3 ) (r£,ef,p£), 

are defined by the forward flow of the following equations: 

ri = ~rieiFi(n,y,ei,pi), (8.18) 

y = (-1 + C ) (?’i))(l + 4>{v)) + (/?(! + Ml) + 0 {ri)){ 1 - <j>{y)), 
h = 2 efFi(ri, y, e 3 , pf), 
pi = eip 1 Fi(n,y,ei,pi), 

(8.19) 

and 

r 3 = r 3 e 3 F 3 (r 3 ,y, e 3 , p 3 ), 

y = (1 + 0{r 3 )){l + </>(y)) + (— /3(1 - p 3 ) + 0{r 3 )){ 1 - <f(y))), 
e 3 = —2e 3 F 3 (r 3 , y, e 3 , p 3 ), 

M 3 = -e 3 p 3 F 3 {r 3 ,il, e 3 ,p 3 ), 
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respectively. These equations are equations (16.1311 and (16.1511 . respectively, written in terms of y rather than 
w , where 

A(ri,y,ei,/ri) = (-1 + 0 {n))(l + <j>(y)) + (a + 0 (ri))(l - f(y)), 

F 3 (r 3 , y, e 3 , p 3 ) = (-1 + 0 (r 3 ))(l + f(y)) + (a + 0 (r 3 ))(l - <j>(y)). 


Lemma 8.12. Consider and let 

F= f 1 -( 1 + </>(y)) +a(l ~<f>(y)) r 

J -r i+m-m -m) 

Then the mappings ^ in 18.1b 1) and take the following form: 

£r : r i = ri(l - Ee i + C>(nei)), 

£± = ei(l + 2Ee\ + 0(riei)), 
A*r = Mi(l + Fe\ + 0(ne i)), 

and 

t 3 '■ = r 3 ( 1 + Ee 3 + 0(r 3 e 3 )), 

e 3 = e 3 (l ^ 2Ee 3 + 0(r 3 e 3 )), 
ht = M3(l - ^e 3 + 0(r 3 e 3 )). 


Proof Consider and (18.181) (the analysis for £ 3 " is identical). Since /3 < 0 we have that y < 0 for 

y € [— 1 , 1 ] for n = yi = 0 and hence we can replace time by y by dividing the equations for f\, and 

pi by y. The point (ri,ei,pi) = 0 is a fixed point of these equations. Solving the second order variational 
equations gives the desired result. □ 

We then have: 

Lemma 8.13. Consider m- The Poincare map P 3 : (r 3 ,e 3 ,p 3 ) 1 —► (rij", £ 3 ", p 3 ) takes the following 
form: 

P 3 ■ = r 3 ( 1 + 2p 3 + A// 2 r 3 + 0(r 3 (e 3 + p 3 + r 3 )), (8.20) 

e 3 = e 3 (l - 2(2 p 3 + A Il2 r 3 ) + 0(r 3 (e 3 + p 3 + r 3 )), 

hj = Ms(l - (2p 3 + Ajj 2 r 3 ) + 0(r 3 (e 3 + p 3 + r 3 )), 


where A jj 2 is given by m- 

Proof. We use Lemma f 8 . 121 and Lemma T 8 .Ill □ 

Remark 8.14. Note that to leading order 18.2(A) is independent of E and hence of f, the regularization 
function. Hence <p does not induce bifurcations in the transition from sufficiently large limit cycles (meaning 
v sufficiently small) in chart K 2 to the 0 ( 1 ) limit cycles. 

Next, we solve for fixed points of P 3 and obtain: 

Proposition 8.15. Suppose that A // 2 in 13.81) is non-zero. Then for p 3 , e 3 sufficiently small and 
p 3 Ajj 2 < 0, the mapping P 3 has a locally unique family of fixed points: 

i ~ 3 = — 2 p 3 A jj 2 + 0(p 3 (e 3 + p 3 )). (8-21) 

The family of fixed point of P 3 corresponds to a C k -smooth family of periodic orbits which are attracting 
(repelling) for Aii 2 < 0 (Au 2 > 0 ;. 

Proof. Suppose that ( r 3 l p 3 ,e 3 ) is a fixed point of P 3 . Then since p = r 3 p 3 and e = r 3 e 3 we solve for 
r 3 = r 3 (p 3 , e 3 ) by setting r£ = r 3 in (18.2011 . This gives 


2 p 3 + A// 2 r 3 + 0 (r 3 (e 3 + p 3 + r 3 )) = 0 . 
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We solve this equation by the implicit function theorem and obtain (18.21[) . The statement about stability 
follows from the fact that the sign of A // 2 determines the sign of d r 3 r^ = r 3 (A // 2 + 0(e 3 + p 3 + r 3 )).D 

The periodic orbits in chart K 2 , described by Proposition 18.31 that are due to the Hopf bifurcation, are 
locally unique since they are obtained by the implicit function theorem for P 2 = x 2 =0. These orbits can 
be continued all the way up to the section Aj, defined in (16.111) (see also Remark l 8 Al) . The periodic orbits 
that are 0(1) with respect to O(e), described by Proposition 18.51 are also locally unique by virtue of the 
implicit function theorem. Therefore by setting £3 = v 2 , corresponding to the section hff (16.lip , and taking 
v sufficiently small, we obtain: 

X 2 = W 1 , y = 1, p 2 = - : ^F r 2 +< 9 ( 7 - 2 ), (8-22) 

using (16.71) and (18.211) . Therefore we can conclude that the periodic orbits described by (18.211) coincide with 
the locally unique ones in chart K 2 described by Proposition 18.31 Similarly, setting r 3 = p, corresponding to 
section A + (14.191) . shows that the periodic orbits in (18.211) coincide with those in Proposition l8.51 where these 
are defined. This gives a (7 fe -smooth and locally unique family of periodic orbits as described by Theorem l 8 . 8 l 
case 7/J. 

Remark 8.16. We believe that the application of the directional charts Ki ,3 in this section to describe 
the Poincare-mapping P e is a novelty. The coordinates of k \ } 3 enabled us to connect the periodic orbits in 
chart H 2 with the larger periodic orbits without the need for careful estimation. This is a general advantage 
of the blowup method and the phase directional charts ^ 1 , 3 . Having said that, it might be possible to prove 
the connection of the limit cycles in k 2 with the larger limit cycles in Proposition 1 8.5\ by working in chart n 2 
alone. To do this one would, however, have to perform a careful estimation of the function D in {8.40 . The 
following remark, Remark \8.17\ contains a discussion of this issue. 

Remark 8.17. By J 6 . 7\) it follows that the family of fixed points of P 3 described in Provosition \8.1 5l 
intersects Aj in 18.22)) . By Theorem I A. (SI case //J, this fixed point of P 3 corresponds to a periodic orbit 
obtained from Provosition \8.3) for a corresponding value of the energy constant h. Hence the value of p 2 in 
18.22)) must agree with the value given in m> - This is a corollary of Theorem, 1 8 .t& case //|. One may 
expect that there could be a more direct way of showing this. We now outline a formal derivation of the result 
by approximation the integrals in D r2 and D ^ 2 in 18.9)) . 

Since X 2 = in 18.221) we take h = To compute D r2 18.6) ) we first substitute dt = ((1 + </>) — 
/3(1 — (f>))x 2 dy and integrate from y (j to yWe then split this integration into (a) an integration from y = y q 
to y = 0 and (b) an integration from y = 0 to y = yjf. We then ignore the contribution from the region 
of regularization and simply set <f> = 4=1 Hi the integrations (a) and (b), respectively. We apply the same 
approximation to the Hamiltonian H(x 2 ,y) and obtain the value ofyQ and yfrom the equation H(0,y) = h. 
Combining this gives the following approximation of D r2 : 

D T2 ~ Zv zp a ( a ^~ + ^ + ) + + X- + a(C + + X + ))) . 

which is valid for v small. Hence 


D r2 


Ajj 2 
J/3 ’ 


from 13.8)) . For D^ 2 , as defined in \ 8 . 7| ), we use 18.81 ) and approximate X 2 by the constant value v 1 . This 
gives 


D 


U 2 


2 

V 


Hence, from m, we have 


92 


Ajj 2 

2v 2 


r 2, 


which agrees with 18.221) . We have not pursued a rigorous result of this kind. 
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8.3.2. Connecting limit cycles for case VI 3 . Krupa and Szmolyan [2T] describe the classical canard 
explosion phenomenon, as observed in the van der Pol system. They prove that the periodic orbits within 
their chart k 2 belong to a smooth and unique family of local periodic orbits that also includes 0(1) periodic 
orbits that arise from their canard explosion. Their proof involves careful estimation on the dependency 
of the function _D(r 2 , y 2 , h), as defined in (18.41) . on the distance to the singular canard (measured by the 
energy constant h). It seems plausible that a similar analysis could be performed to our system. However, 
the situation here is complicated by the fact that our Hamiltonian function H depends non-trivially on 
the regularization function </>. The statements in Theorem 18.81 case VTJ, therefore summarize the previous 
results in Proposition 18.31 and Proposition 18.61 Instead we conjecture on the connection of small periodic 
orbits in Proposition 18.31 with the larger ones in Proposition 18. 61 as follows: 

Conjecture 1. The two families of periodic orbits in V/J belong to the same family of locally unique 
periodic orbits. 

The limit cycles in Proposition 18.31 and Proposition 18.61 do not seem to be present in the PWS case VI 3 . 

8.4. Saddle-node bifurcation. We conclude this section with the following main result 

Theorem 8.18. Suppose j 7 .1 | ) and \5.ff . Then for e sufficiently small: 

J/|: There exists an open set of regularization functions such that the periodic orbits in Theorem Iff. <81 
case 112 , undergo at least one saddle-node bifurcation. 

V /J: Suppose, in addition, that Conjecture [7] holds. Then there exists an open set of regularization func¬ 
tions such that the periodic orbits in Theorem. 1 8.81 case H/J, undergo at least one saddle-node 
bifurcation. 

Proof. A corollary of Theorem 17.61 is that we can always achieve a 2 A // 2 < 0 or a 2 Av "/ 3 < 0- The result 
therefore follows from the statements in Theorem 18.81 and Conjecture [I| □ 

9. Numerics. In this section we illustrate the results in Theorem 18.181 and provide further support for 
Conjecture [l] by computing limit cycles for two model systems, for cases //| and VJJ. 

9.1. Case ITj. In this section we consider the regularization of the following model system for case II 2 ' 

X+{X "V) = {~x+~2 7 J) ’ *" (W) = ((x-J-^x-p) 2 ) ’ ^ 

corresponding to the following parameters: 

5 = -l,a = l,0= -1, C + = -7, C = - 6 , V ± = ±2, X ± = 0, fi = 2, 
in (12.121) and (12.131) . The constant Ajj 2 given in (13.81) takes the value: 


A n 2 = -6. 


According to Proposition 13.61 the limit cycles of the PWS system are therefore all stable. According to 
Theorem 18.81 the 0(1) (with respect to e) limit cycles of the regularized system are also stable. 

We consider two regularization function^]: 

4 >\y) = v for ye (-1,1), (9.2) 


<t> c {y) = 7jy-\tf for ye (-1,1). 

From (17.131) we obtain j/q = 0 in both cases. Inserting the corresponding values of 
(j>"( 0) = 0, (p 3 t H = ^ 3 )(0) into (17.141) and (17.151) gives the following values for and a 2 : 

a 4 , h = 13?' 2 + Oirf), 

, 1 

a 2 2 ’ 

26 9 

y C 2,H = y r 2 + Ofrl) = (8 + -)r 2 + 0(r 2 2 ), 



(9.3) 

^( 0 ), 4>2,h = 

(9.4) 


7 The superscripts l and c refer to linear and cubic, respectively. 
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Figure 9.1. Periodic orbit amplitudes as a function of the parameter [12 for r 2 = 0.01. The full line shows the result of 
using the linear regularization function El) , the dotted line shows the result of using the cubic regularization function El- 
As can be seen, the linear regularization function induces a saddle-node (SN) bifurcation, in agreement with Theorem \ 8 . 18 l 


Since a l 2 Ajj 2 < 0 we have from Theorem 18.181 that the linear regularization function (19.21) introduces a 
saddle-node bifurcation. We demonstrate this as follows. 

Using the numerical bifurcation software AUTO we continued the periodic orbits in the two regulariza¬ 
tions of II 2 ■ In Fig. 19.11 we show the amplitude (measured as max(y)) of the periodic orbits as a function 
of the parameter /x 2 for r 2 = \fe = 0.01. The full line shows the result of using <j>i , as given in (19.21) . while 
the dotted line shows the result of using </> c , as given in (19.81) . The linear regularization function introduces 
a saddle-node bifurcation, in agreement with Theorem 18.181 On the other hand, the cubic regularization 
function does not introduce any saddle-node bifurcation. 

The Hopf bifurcations were numerically found to occur at 

ft l 2 H « 0.1298, fi c 2}H « 0.0866, 

which are in good agreement with (19.41) for r 2 = 0.01. Also, in agreement with Proposition [53J we observed 
that the two family of limit cycles agree for larger values of /. x 2 since the limit cycles for both regularizations 
must be 0(e)-close to the limit cycles of the PWS system for p = 0(1). 

Remark 9.1. We obtained system m by fixing m and m and the parameters 8 , a, (3 and solving 
a 2 = —a 2 for and r 7 ± setting = 0 for simplicity. 

9.2. Case VJJ. In this section we consider the following model system for the case VI 3 : 

X+ = C-|I)' X_ = (-2(z-rfl (*-„)*)' <“> 

Here 

8 = 1, a = -1, /3 = 2, C + = C” = 0, v + = 0 , ry - = 1 , x ± = 0 , U = 1 . 

We sketch the PWS system in Fig. 19.21 for pt = 0. It is very similar to Fig. 12 in [21]. However, as opposed 
to HU we have included the cubic term in A' + which gives rise to an invisible tangency at (a;, y ) = (1, 0) and 
a return mechanism from £"*" to £“. 
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Figure 9.2. Sketch of the PWS system 19.51) for p = 0. The triple-headed arrows within E indicate the direction of the 
sliding vector field. 


We then regularize X in (19.51) using the cubic function in (19.51) . Since H = 1 > 0, we can apply 
Proposition 17.31 to the regularized system and conclude that the system has an equilibrium (17.111) which 
undergoes a Hopf bifurcation at 


M2 = H c 2 ,h = ~7^ r 2 + °( r 2)- (9-6) 

The first Lyapunov coefficient is obtained from (17.151) : 

« = ~^ 2 +0(rl). (9.7) 

Since a 2 = — ^ < 0 we conclude that the periodic orbits are attracting for r 2 (and hence e) sufficiently small 
and appear for /x 2 > H . 

This example also has a maximal canard (see Theorem l6.4D . The parameter value at which it occurs is, 
from (16.221) . given by 


A* 2 ,c = Try r 2 + 0{rl) « —0.07806r 2 + 0(r|), (9.8) 

90i, c 

In the last expression we have used (16.101) to obtain 

« 1.4233. 

There is a related example for the case VI 3 in [22] on p. 2169 (after reversing time and reflecting x 1 —> — x) 
with the same values of 8 , a and /3. The system in [22] has = q=l as the only non-zero coefficients in 
(12.111) . The reason for modifying the system given in [22] is that their system gives a 2 = 0 from (17.151) , for all 
regularization functions <f>. In fact a detailed calculation shows that a = 0. The example in [22] is therefore 
codimension two for the regularization. 

In Fig. 19.31 we have used the numerical bifurcation software AUTO to track the amplitudes of the limit 
cycles of (19.51) emanating from the equilibrium (17.21) . We considered r 2 = 0.1. The amplitude of the limit 
cycles is now measured in Fig. 19.31 using max(s) instead of max(y) used above. This proved to be more 
illustrative in this case. A dramatic increase in amplitude is seen near /i 2 « —7.836574 x 10 -3 . In Fig. 19.41 
we have illustrated three different limit cycles within the original (a;, y)-plane. The largest limit cycle looks 
like a canard. The three limit cycles occur for the following parameters: 

= —7.85 x 10" 3 , 

= -7.8365738 x 10" 3 , 

= -7.8365737 x 10" 3 . 
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Figure 9.3. The amplitude of the limit cycles, measured in terms of max(x), as a function of p 2 - The amplitude explodes 
due to the presence of a maximal canard around p 2 ~ —7.8366 x 10 —3 . The flat part beyond the canard explosion gives rise to 
the canard-like limit cycles similar the one shown in Fig. \ 9.f\ 


The difference between the last two parameters is 10 11 . The dramatic increase of amplitude is due to the 
canard explosion phenomenon described in Pronosition l8.6l Numerically we found the following canard value 


l4 tC = -7.8365738 x 10“ 3 . 


(9.9) 


This value is in good agreement with (19.81) for r 2 = 0.1. Note that in comparison to the classical canard 
relaxation oscillation in the van der Pol system, the duck’s head and chest are in our case, not due to motion 
along a curved slow manifold. Instead they are due to regular motion within £± following the regular vector 
fields A , respectively. It is the motion along the slow manifold that creates the straight back of the duck. 
Also in the present case these different types of motions occur on an identical time-scale. There is a slow-fast 
behaviour but it is hidden and only visible through the scaling y = y/e. 

Now we replace the cubic regularization function in (El by the following septic C 1 regularization 
functiorU: 

&% y ) = + T ^ ii 5 + y , ye (-1,1). (9.10) 

This regularization function has been constructed so that < 22 , using ( 17 . 151 ) . becomes 



This value is just the negative of the previous value in E3- Hence periodic orbits emanating from the 
Hopf bifurcation are repelling and appear for y < /j,2 ,h where now 

= ^ 2 ,h — ~ g\/e + ^( e )- (9-11) 

The canard value y 2 , c a l so changes and becomes 

/4 iC « —0.12188^ + 0(6). 


8 The superscript s now stands for septic. 


(9.12) 
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Figure 9.4. Three different limit cycles. Due to the canard explosion phenomenon (see Fig. \M> 
parameter between the largest and the second largest limit cycle is extremely small: 10“ 11 . 


the difference in 


We again use AUTO with r 2 = 0.1 to continue periodic orbits from the Hopf bifurcation at /z = fi s 2 H (19.111) . 
We obtain the bifurcation diagram in Fig. 19.51 As opposed to Fig. 19.31 we now observe a saddle-node (SN) 
bifurcation, which occurs before the canard explosion phenomenon. 

Remark 9.2. Fixing the values of S, a and fd it is straightforward to construct a family of model 
systems for VI 3 where the Lyapunov coefficients corresponding to regularization functions <jr and <j> c (see 
H9.2\) and Iff.ffl ) ) have opposite signs: < 0. A simple example is the following: 


X + (x,y) = 


l + 7x 
t + 8x 2 


X (a;,y,/z) = 


— 1 + 6a; 

— §(a; - /z) + 8(a; - nf 


where 


a 2 



a 9 = 


7 

32' 


However, we have not presented the details of this case since the saddle-node bifurcation occurs very close to 
the canard value and therefore it is not as clearly visible as the saddle-node in Fig. Iff. 51 

Remark 9.3. This “duck” part in Fig. \9.f\ is not covered by our results. However, it seems very plausible 
that the results in [2 if can be extended to this case too. 

10. Discussion and Conclusions. In this paper, we have considered the regularization of the codi¬ 
mension one two-fold bifurcation in planar PWS systems. The PWS two-fold bifurcation is dynamically 
very interesting as it may include singular canards, pseudo-equilibria and limit cycles. Using the blowup 
method of Krupa and Szmolyan m , we continued these objects into the regularization and we related the 
PWS bifurcations to standard smooth bifurcations. Perhaps most interestingly, we were able to show that 
the regularization can induce saddle-node bifurcations of the limit cycles. The results were illustrated by 
numerical examples. 

There are two questions that emerge from this work that we feel are worthy of further discussion. What 
light can regularization shed on the original PWS system? Is the introduction of saddle-node bifurcations a 
necessary consequence of regularization? 

For the first question, it is clear that singular canards and pseudo-equilibria of the PWS system are 
limits of equivalent objects in the regularized system. Similarly, limit cycles in the PWS two-fold case II 2 
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V 2 


Figure 9.5. This diagram shows the amplitude of the limit cycles, measured in terms of max{x), as a function of jjrj 
using the septic regularization function 19.10\) to regularize the system S 9. 51 ). The amplitude explodes due to the presence of a 
maximal canard around p,2 ~ - 1.22369 X 10 2 . 


are limits as e —> 0 of limit cycles in the regularized case //J, at least “macroscopically”; the saddle-node 
bifurcations occur “microscopically” within chart K 2 - In comparison, the PWS case VI 3 is more singular. 
It possesses backwards and forwards non-uniqueness of orbits due to the presence of stable and unstable 
sliding. In particular, it is possible to identify closed “singular cycles”, reminiscent of singular cycles in slow- 
fast systems such as the van der Pol system (see Fig. m- Our analysis showed that these singular cycles are 
limits of periodic orbits of the regularization. Interestingly, the quantity A yi 3 defined in (18.111) only depends 
on the PWS system, giving us an insight into the stability of a very singular object. The limit cycles of the 
regularization undergo a canard explosion phenomenon which gives rise to a very rapid amplitude increase 
of local periodic orbits. This can lead to global limit cyles as it was shown in section F!~?l and Fig. 19.41 

The second question is much broader. In this paper, we have considered planar two-folds, subject to the 
Sotomayor and Teixeira [30] regularization. We have shown that the criticality of Hopf bifurcations depends 
on the regularization function and generically it is possible to induce saddle-node bifurcations by varying 
the regularization function. But we have not shown how many saddle-node bifurcations may exist. Perhaps 
there are other PWS systems where the regularization does not induce bifurcations. Or there may be systems 
where other types of behaviour occur upon regularization. In addition there are other regularizations that 
could be considered. 
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